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ABSTRACT 


Search and Detection Theory is the overarching field of study that covers many scenarios. 
These range from simple search and rescue acts to prosecuting aerial/surface/submersible 
targets on mission. This research looks at varying the known discrete and Bayesian algorithm 
parameters to analyze the optimization. It also expands on previous research of two searchers 
with search radii coupled to their speed, executing three search patterns: inline spiral search, 
inline ladder search, and a multipath ladder search. Analysis reveals that the Bayesian search 
and discrete search work similarly, but the Bayesian search algorithm provides a more useful 
output in location probability. Results from the continuous search were similar to previous 
research, but variance in time to detection became more complex than basic increasing or 
decreasing ranges. 
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CHAPTER 1: 

Introduction 


The process of finding the location of one or more unknown subjects is known as Search and 
Detection Theory (SDT) and is of great interest to U.S. Navy (USN). While there are many 
applications, both military and non-military, SDT is of particular import to Anti-Submarine 
Warfare (ASW). This is due to the desire to locate adversarial submarines because the nature 
of submarines is to use their covert operating patterns to conduct Intelligence, Surveillance 
and Reconnaissance (ISR), Anti-Surface Warfare (ASUW), Strike Warfare (SW) and Special 
Operation Forces (SOF) missions. 

With the advent of high end diesel-electric submarines, increased tensions and activity 
from capable nuclear submarines (Russia and China), and the proliferation of submarine 
technologies to previously non-capable navies, the need for better approaches to SDT is 
evident. The proliferation of these submarine technologies, and the rise of autonomous 
vehicles such as unmanned underwater vehicles (UUV) and unmanned surface vehicles 
(USV), both as preprogrammed searchers and potential targets, dictate that the parameters 
of SDT algorithms must be investigated. 


1.1 History 

The Operations Evaluation Group (now known as the Center for Naval Analyses (CNA)) [1] 
was created during World War II (WWII) to enlist the help of mathematicians and scientists 
to assist in developing new ASW tactics and strategies. This small group of MIT scientists 
worked with ASW Operations Research Group (ASWORG) directly in the field for optimum 
observations and real-time solutions [2]. These scientists used the burgeoning search theory 
to improve U-boat search and neutralization, especially with favorable patterns and operating 
characteristics for patrol aircraft. 

After the war, search theory was extended to other naval operations such as searching 
for the USS Scorpion (SSN-589), lost ordnance, and tracking Russian ballistic missile 
submarines [3]. The United States Coast Guard (USCG) developed a computer-aided 
search program (CASP) [4] for missing boats and mariners lost at sea. 


1 




Dr. Bernard Koopman of the WWII ASWORG developed the historical basis for search 
theory by dictating the following basic principles. 

1. Quantify the probability of the sensor capability to detect the target. 

2. Classify the movement of the searchers and target. 

3. Analyze the constraints on the searcher’s forces and resources 

4. Develop the highest probability of detection with the greatest economy of forces [5]. 

1.2 Applications and Previous Research 

The applications for SDT are numerous, so that while the context of this research is limited 
to an ASW submersible search problem, the basics can be extended to various problem sets 
and models. Research conducted at Naval Postgraduate School (NPS) has focused on naval, 
non-naval, and other Joint Warfare (JW) applications. Research has also been conducted at 
non-military institutions for civilian applications. 

1.2.1 Some Previous Research 

1. Zachary Lukens’s 2016 thesis [6] investigated SDT applied to a single target in the 
South China Sea and two searchers using cookie cutter sensors utilizing inline spiral, 
inline ladder, and multipath ladder searches. His research found that increasing either 
the detection radii or speeds decreased the time to detection and that spiral search 
was faster than ladder searches at large radii, but both ladder searches were faster than 
spiral search for small radii. 

2. Volkan Sozen’s 2014 thesis [7] used SDT applied to Unmanned Aerial Vehicles (UAV) 
patrolling a border. He modeled the search as one UAV searching the area, two UAVs 
searching the same area, and two UAVs independently searching two equally split 
sub-areas. His research found that detection increased with the second searcher 
and splitting the searchers into independent sub-areas, the border was covered more 
efficiently. 

3. Dr. Michael Atkinson [8] (together with Moshe Kress and Rutger-Jan Lange) re¬ 
searched searching for a target discrete location based on receiving of intelligence 
of varying reliability and the time cost incurred based on searching a certain area. 
His research became the problem of looking at the probability of detection based on 
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amount of intelligence gathered weighed against the target escaping due to a delayed 
search. This has direct applications to military and non-military problems based 
on the information quality. The study concluded that high intelligence credibility 
had better results in waiting to attack, while low credibility had better results for 
immediate attack. 

4. H. Zhou and H. Wang in [9] examined optimum paths for a solitary searcher looking 
for a solitary target using spiral, ladder, loop and square spiral paths.They showed 
that a stationary searcher is optimally placed in the center of the search area. It also 
showed that, with a moving searcher, the ladder pattern gives the most likely detection. 
In [10], they examined optimum searcher placement seeking a target travelling back 
and forth between a detectable operating area and an undetectable hiding area. They 
showed closed form solutions for determining mean time to detection based on the 
placement of the searcher in a given path from hiding to operating area or within 
the operating area itself. Their results showed how to distribute sensors ("searchers") 
among the routes and the operating area. They extended their study to include sensors 
in detecting a target moving between a hiding area and an operating area [11]. Based 
on the average time to detection, they were able to evaluate the performance of placing 
the sensor(s) to monitor various travel pathway(s) or to scan the operating area. 

1.2.2 Applications 

The USN has utilized SDT for 

• ASW 

• ASUW 

• Anti-Air Warfare (AAW) 

• SW 

• Mine Warfare (MIW) 

• other operations. 

Other branches and JW applications include 

• ISR 

• Counter-drug operations 

• Counter-terrorism operations 
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• Ballistic Missile Defense 

• SW and Artillery Fires operations 

• Seareh and Reseue. 

Multiple eivilian applieations are also prevalent in researeh and development, to inelude: 

• Search and Rescue 

• Quality Control 

• Resource exploration and seareh (oil, minerals, fishing, ete.) 

• Mieroeeonomies (resouree to market matehing). 
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CHAPTER 2: 
Discrete Search 


2.1 Theory 

The following sections outline the basic concepts and theory behind the discrete search 
theory. 

2.1.1 Terms 

A - Total search area (no searches conducted outside of A) 
n - Total number of equal area subdivisions (boxes) of A : (i £ {l..n}) 

Pi - /'[tgt in /th box] probability of target present in zth box 

qt - /’[miss detection |tgt in /th box] probability of failing to detect target in /th box 
given it is present in /th box 

ki - number of glimpses (discrete time interval) with search was conducted in /th box 

n 

m - total number of glimpses: m = Yj 

i=\ 

n 

By the concept of total probability, Y Pi = 1 • 


2.1.2 General Concepts 

The basic assumptions for the discrete search model are [12]: 

• Single stationary target 

• Single searcher with imperfect sensor (qi) 

• Total search time limited (m glimpses) 

• No time lost transitioning from one box to another 

• Each search conducted independent of other glimpses 

• No false positives 

• Non-adaptive search. 

The probability of the target not being detected in the /th box on a given glimpse is 
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P[miss detection] = P[miss detection | tgt in /th box]P[tgt in /th box] = = piqi and 

the detection probability on a given glimpse for the ith box is /’[detection] = /’[detection | 
tgt in /th box]/’[tgt in /th box] = = /’/(I - qd- After ki glimpses, the probabilities of 

missing or detecting in the /th box are /’j^‘ = /?/^f‘ and = p/(l - qf‘), respectively. This 

leads to [13], by the total probability concept, the likelihood of total failure and detection as 

n 

PM = '^Piqf‘ ( 2 . 1 ) 

i=\ 

and 

n 

PD = YjPi^^-^h ( 2 . 2 ) 

i=\ 


2.1.3 Model Overview 

The greedy algorithm is a myopic search model that is designed to attempt to optimize 
Equation 2.1 in order to streamline the detection based on a limited number of resources 
(time, searchers, etc.). For each glimpse, the probability of detection is calculated given 
previous glimpses have failed to detect the target. This change in detection probability [13] 

A = (1 - - (1 - qf)pi = (1 - qi)q\'pi (2.3) 

is used to select the maximum detection probability and the search is conducted in that box, 
increasing k,- by 1 for that box. This updates the probability for that box alone. If the target 
is not found, the search proceeds to the next glimpse using the same formula. This will be 

n 

repeated until the target is found or 2 h = m. The cumulative detection probability is the 

1=1 

sum of the maximum A for at every step or Equation 2.3. 

2.1.4 Example 

The example in Table 2.1 from [12], outlines the "greedy" search algorithm. The target 
is in one of three boxes, with associated location and conditional miss probabilities. The 
probability of detection for each cell is calculated based on finding the target on the first 
glimpse with no previous looks (k = (0,0,0)) so that q'^' = 1. The highest probability is 
in the first cell making the qf' < 1 reducing the probability on the next round of glimpses. 
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The cumulative probability after six glimpses is the sum of the maximum probabilities at 
each round, Pd = .25 + .125 + .1 + .08 + .064 + .0625 = .6815. 

Table 2.1. Optimal Sequence of Discrete Looks. Adapted from [12]. 


{ki,k2,k2,) 

/>i(l - 

P2{1 - q2)q\^ 

/>3(1 - 

(0,0,0) 

.25 

.08 

.1 

(1,0,0) 

.125 

.08 

.1 

(2,0,0) 

.0625 

.08 

.1 

(2,0,1) 

.0625 

.08 

0 

(2,1,1) 

.0625 

.064 

0 

(2,2,1) 

.0625 

.0512 

0 


P\ 

P2 

P3 

.5 

.4 

.1 

q\ 

q2 

qs 

.5 

.8 

0 


2.2 Problem and Model 

This basic model for discrete search is based solely on the conditional miss probabilities 
{qt) and the location probabilities (p/). These parameters will be varied and results shown. 

2.2.1 Parameter 

The location probabilities will be varied according to Figure 2.1 and the detection probabil¬ 
ities will be varied according to Figure 2.2. The uniform location and conditional detection 
probabilities are to set up the base case to compare other distribution combinations too. 

The base location probability per box, = ^ = as shown in the uniform case. Two 
linear stratified location probability densities are used to simulate the target being placed 
preferentially along some line of longitude (small scale) for the first distribution or some 
line of constant longitude/latitude line. Each of these increase by a factor of 10 over the 
search area in the direction of greatest change. Three Gaussian distributions are used to 
simulate the likelihood of the target having a well-localized location probability prior. They 

( x - xq )^ (y-yo4 

follow a distribution function p{x, y) = —\=e . The centered Gaussian has 

cryln 

different criat and cr/o„^. The other two Gaussian plots are centered on a corner and side 

with CTiat = (Tlong- 
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The conditional detection probability densities are uniform as a base sensor that works 
the same in all locations qi = .3, giving a 70% of detection at each location. The two 
stratified detection probability densities are linear stratified along latitude or longitude from 
qi = (.063, .63). This simulates the sensor performance based on some geographic profile 
(i.e., change in bottom type for SONAR). 


Figure 2.1. Location Probability Distributions pi 



From top left clockwise: Uniform dist., Linear stratified along longitude, Linear 
stratified along Long/Lat, Gaussian dist., Gaussian centered in corner, Gaussian 
centered on side 


Figure 2.2. Conditional Miss Probability Distributions qi 



along latitude 
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2.3 Results 

The following sections show the results for six different location probability densities and 
the three conditional search probability densities. Each combination was examined for 
n = 50 glimpses (half the number of available boxes). 

Since there are 18 possible combinations each with various tables and figures, the following 
will show a small sample of the results. All of the figures and tables will be presented in 
Appendices A and B, respectively. 

2.3.1 Linear Stratified Location 

Figures 2.3, 2.4, and 2.5 show the remaining location probability density for the three 
conditional search probability densities after 50 glimpses. Also shown are the boxes that 
were selected to search by the greedy algorithm. 

For the uniform conditional search probability in Figure 2.3, unsurprisingly, the selected 
boxes start with the highest pi and slowly work to the left as p,- lowers. 


Figure 2.3. Stratified Location and Uniform Search Density Results 
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Remaining probability of detection and boxes search based on stratified p,- and 
uniform qi 


For the longitudinal stratified conditional search probability in Figure 2.4, the selected 
boxes also started with highest pi and slowly work to the left as pi lowers. This is slightly 
unexpected since a better probability of detection exists in areas of lower target location 
probability. But the factor of 1 - qi is 93.7% in the western side and 37% in the east. This 
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coupled with the pi being a factor of 10 higher in the east leads it to be a higher probability 
of detection even though the sensor performs worse there. 


Figure 2.4. Stratified Location and Longitudinal Stratified Search Results 



Remaining probability of detection and boxes search based on stratified p,- 
and linear stratified qi 


Figure 2.5. Stratified Location and Corner Stratified Search Density Results 



Remaining probability of detection and boxes search based on stratified pi and 
corner stratified qi 


For the corner stratified eonditional seareh probability in Figure 2.5, the seleeted boxes start 
in the southeast due to the juxtaposition of high pi and low qi. The path moves mostly 
along diseretized radial ares eentered on the southeast eorner. Whieh box is ehosen via the 
greedy algorithm is dependent on when 1 - qi will push a probability of deteetion higher 
for a lower p,-. 
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For the longitudinal stratified pi location probability density, the following results were 
found based on running the search for up to 500 glimpses. A data set using a step size of 
5 is shown in Figure 2.6 along with associated fit lines as calculated from MATLAB’s fit 
algorithm based on: 

fix) = 1 - e-‘’f 

Each of the fit lines are color coded to the data set with black for uniform, red for longitudinal 
stratification, and green for corner stratification. As the number of glimpses increases, the 
probability approaches 1, with the uniform case increasing faster than the other two. 

Figure 2.6. Stratified Location Cumulative Probability and Fit Lines 

Pd vs Number of Olimpses and Exponential Fitted Curves 

1 
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Cumulative probabilities of detection based on the number of glimpses. In¬ 
cluded are the fit best fit line 
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Table 2.2. Parameter Variation for p, for Longitudinal Stratification 


Parameter 

b 

r2 

SSE 

RMSE 

P strut j^uni form 

P strat->Qlong 

P strut jQcorner 

0.013208 

0.0089597 

0.011187 

0.99908 

0.99920 

0.99808 

0.00065853 

0.00066179 

0.0013515 

0.008554 

0.0085751 

0.012254 


Best fit line parameters for a longitudinally stratified location probability density 
of the form f{x) = 1 - 


2.3.2 Gaussian Corner Location 

The remaining loeation probability densities for the three eonditional seareh probability 
densities after 50 glimpses are shown in Figures 2.7, 2.8, and 2.9. Also shown are the boxes 
that were seleeted to seareh by the greedy algorithm. 


For the uniform eonditional seareh probability in Figure 2.7, unsurprisingly, the seleeted 
boxes start with highest p, and slowly work outward as p/ depletes. A seeond look begins 
oeeurring in those higher probability boxes prior to searehing all boxes. 


Figure 2.7. Gaussian Corner Location and Uniform Search Density Results 



Remaining probability of detection and boxes search based on Gaussian pi and 
uniform qi 


For the longitudinal stratified eonditional seareh probability in Figure 2.8, the seleeted 
boxes also started with highest p,- and slowly works outward as p,- lowers. This is slightly 
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unexpected since there is better probability of detection in some areas of lower target 
location probability. As the probability of detection decreases on subsequent looks, the 
lower qi increases the likelihood of finding the target to the west of the initial location 
locus. From the Remainder Probability graph in Figure 2.8, the higher location probability 
remains dispersed in the northwest corner albeit non-uniformly. 


Figure 2.8. Gaussian Corner and Longitudinal Stratified Search Density Re¬ 
sults 
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Remaining probability of detection and boxes search based on stratified pi and 
linear stratified qi 


For the corner stratified conditional search probability in Figure 2.9, the selected boxes start 
in the northeast due to the high p, overriding the high qi. The Remainder Probability graph 
in Figure 2.9 shows some areas 2-3 times greater than other areas, especially around the 
original locus. Some areas will be searched multiple times before every box is searched 
once. 
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Figure 2.9. Gaussian Corner and Corner Stratified Search Density Results 



Remaining probability of detection and boxes search based on stratified pi and 
corner stratified qi 


For the Gaussian Corner pi location probability density, the following results were found 
based on running the search for up to 500 glimpses. A data set using a step size of 5 
is plotted in Figure 2.10 along with associated fit lines as calculated from MATLAB’s fit 
algorithm based on: 

m = 1 - e->’\ 

Each of the fit lines are color coded to the data set with black for uniform, red for longitudinal 
stratification, and green for corner stratification. As the number of glimpses increase, the 
probability approaches 1, with the uniform case increasing faster than the other two. 
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Figure 2.10. Gaussian Corner Cumulative Probability and Fit Lines 



Cumulative probabilities of detection based on the number of glimpses. In¬ 
cluded are the fit best fit lines 


Table 2.3. Parameter Variation for pi for Corner Gaussian Location 


Parameter 

b 


SSE 

RMSE 

PGC ^Quniform 
PGCi^long 
PGCiQcorner 

0.03117 

0.01869 

0.01441 

0.99683 

0.9952 

0.99523 

0.0013941 

0.0027014 

0.0029724 

0.012446 

0.017325 

0.018173 


Best fit line parameters for a Gaussian Corner location probability density of the 
form f{x) = 1 - 


2.3.3 Cumulative Detection Probability 

The total probability of detection for the after 50 glimpses is shown in Table 2.4. The 
Gaussian p,’s have higher probabilities, which is not unexpected since they have a much 
higher concentration on a few points. This allows for multiple searches over a smaller 
number of points instead of 50 different boxes. The other result seems to be that higher 
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probability results when pi matehes up to 1 - qi as it does when pstratified has a higher 
probability of deteetion when paired with qcomerstrat instead of qiang^trat. 


Table 2.4. Po after 50 glimpses 


Pi 

Uniform Long. Stratified 

Corner Stratified 

Uniform 

0.35 

0.4055 

0.44906 

Long. Stratified 

0.50909 

0.34927 

0.46119 

Corner Stratified 

0.57157 

0.43958 

0.40031 

Gaussian 

0.55932 

0.52965 

0.5919 

Gaussian Side 

0.71881 

0.55738 

0.65176 

Gaussian Corner 

0.79975 

0.63655 

0.54999 


Total probability of detection after 50 glimpses based on various combinations 
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CHAPTER 3: 
Bayesian Search 


3.1 Theory 

The following sections outline the basic concepts and theory behind the Bayesian approach 
discrete search theory. 

3.1.1 Terms 

The same terms from Chapter 2 apply along with the following random variables: 

Li - Target is located in box i 
Mi - Target is not detected in box i 

/?o - P[tgt not in A] probability of target not in search area 

n 

By the concept of total probability, Y, Pi = Pa and \ - Pa = Pq- 

i=\ 


3.1.2 General Concepts 

We use the same basic assumptions of Section 2.1.2 with an added assumption that the 
target may not be in the search area. The basics of Bayesian search theory is to update the 
probability of the target being located in the fth (p, ) box based on previous failures to detect. 
Searching the zth box and failing to find it will decrease its likelihood of being located there 
and increase the probability of its being located in another box. This posterior probability 
comes from Bayes’ Theorem. 


P{Li I Mi) = 


P{Mi I Li)P{Li) 
P{Mi) 


(3.1) 


The probability of not detecting the target in the z'th box (P(M/)) is currently the only 
uncalculated quantity. It is, by total probability, the sum of the conditional detection 
failures based on the location probabilities. P{Mi) = P{Mi \ Li)P(Li) + P{Mi \ -<Li)P{-<Li). 
P{Li) and -'P(Li) are pi and 1 -pi, respectively. P(Mi \ Li) = qi and P{Mi \ ->Li) = 1 based 
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on the assumption of now false positives. This updates Equation 3.1 to 


p' = P(Li I Mi) = 


qiPi 

(1 - Pi)+Piqi' 


(3.2) 


The probability that the target is located in jth box {j i) must also be updated since it is 
more likely that it is not located in the ith box, 


P(Li I Mi) 


P{Mi I Lj)P(Lj) 
P(Mi) 


(3.3) 


where /’(L,) = pj and E(M/) = (1 - pi) + piqi as found above. P{Mi \ Lj) is the conditional 
probability of failing to detect the target in box i given it is in box j. Since there are no false 
positives, the searcher will always fail to detect when searching in a box where the target is 
not located, making P{Mi \ Lj) = \. This updates the location probability by 




P(Li I Mi) 


Pj 


(1 - Pi) + Piqi' 


(3.4) 


This updates all pi^j including pQ. 


3.1.3 Model Overview 

In the Bayesian search model, the conditional probabilities of detection failure only factor 
into updating the new location probabilities. As such, the search algorithm searches the 
box with the highest likelihood of containing the target and searches. Upon failure to find, 
it updates all location probabilities and repeats until the target is found or the number of 
glimpses has been reached. 

This still follows the Bayesian search algorithm from Chapter 2 in that at each step the 
probability to detect Pdi = ^/(l - qt) is calculated for all boxes. The search is conducted in 
the box with the highest detection probability. After a failure to detect, the p/’s are updated 
via Equation 3.2 and Equation 3.4 [14]. The Bayesian search algorithm then treats the next 
glimpse as the first with the updated location probabilities [14]. 

If po = 0, this will provide the same output as the Discrete Search Method. Eor example, 
the updated location probabilities are just A,- from the Discrete Method divided by 1 - ^/’s 
and then normalized so that the sum is equal to 1. Multiplying 1 - ^,’s back in after 
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normalization will make the same evaluation as in Table 2.1 [14]. 


3.1.4 Examples from Other Research 

The Bayesian method has been used for several real-world searehes. Famously used in USS 
Scorpion, a submarine lost in 1968, it was also used in MV Derbyshire, the SS Central 
America, the lost hydrogen bomb from a B-52 erash, Air Franee flight 447, and most 
reeently flight MH370. Figure 3.1 is the loeation probability density originally eonstrueted 
for the USS Scorpion [15]. 


Figure 3.1. USS Scorpion Location Probability. Source: [15]. 
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The basie steps outlined in previous examples are: 
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1. Construct plausible conjectures as to the cause of the loss of the object. 

2. For each conjecture, construct a location probability density map (F(L,) = p,). 

3. Build the conditional detection probability maps for every location covered under the 
location probability map (F(M/ | L/) = qi). This will be a function of sensor’s ability 
to detect in the given environment and must be built individually based on current 
conditions and sensor dynamics. 

4. Combine the maps in Step 2 and 3 to create an overall detection probability density 
map (p/(l - qi)). 

5. Develop a search path that scans the highest to lowest probability areas. 

6. Revise the location probabilities continuously during the search. If multiple conjec¬ 
tures are being tested simultaneously, each should be updated independently and if 
the search fails to corroborate aspects of a particular conjecture the likelihood of that 
scenario will need to be downgraded [15]. 

One last aspect of the Bayesian search, that is not being explored in this research, is the 
concept of the cost of a search. If, for example, searching one particular box takes longer 
time, requires more resources, or is more dangerous, then that should be taken into account 
as a weight on the search. If the cost function can be estimated as c/, then the new detection 
probabilities come out as [16]. 

3.2 Problem and Model 

This basic model for Bayesian search is based solely on the conditional miss probabilities 
iqi) and the initial location probabilities (poi)- These parameters will be varied and results 
shown. 

3.2.1 Parameters 

The location and conditional miss will have the same variations as presented in Chapter 2 
(Figures 2.1 and 2.2) and they are repeated in Figures 3.2 and 3.3 for the reader’s convenience 
with po = .1 and the detection probabilities will be varied according to Figure 2.2, which 
are shown again here for clarity. 
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Figure 3.2. Bayesian Initial Location Probabilities poi 



From top left clockwise: Uniform dist., Linear stratified along longitude, Linear 
stratified along Long/Lat, Gaussian dist., Gaussian centered in corner, Gaussian 
centered on side 


Figure 3.3. Bayesian Miss Probabilities qt 
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From left to right: Uniform dist.. Linear stratified along longitude. Linear stratified 
along latitude 


3.3 Results 

The following sections show the results for six different location probability densities and 
the three conditional search probability densities. Each combination was examined for 
n = 50 glimpses (half the number of available boxes). 

The tables and figures showing all 18 possible combinations are presented in Appendix A 
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and Appendix B, respectively. 


3.3.1 Linear Stratified Location 

Figures 3.4, 3.5, and 3.6 show the remaining location probability densities for the three 
conditional search probability densities after 50 glimpses. Also shown are the boxes that 
were selected to search by the Bayesian search algorithm. 


For the uniform conditional search probability in Figure 3.4, unsurprisingly, the selected 
search boxes start with highest pi and slowly work to the left as pi lowers. This follows 
the same pattern as the discrete search from Chapter 2. The exception being that the 
upper-leftmost graph of Figure 3.4 now shows the updated location probabilities in Bayes’ 
algorithm instead of merely the next search step as in the discrete search method. This 
follows as the Bayesian method and Discrete method calculate roughly the same probability 
of detection with the exception that the Bayesian search modifies pq = .1 ^ 0.1846 . 


Figure 3.4. Stratified Location and Uniform Search Density Results 



Remaining probability of detection and boxes search based on stratified pi and 
uniform qi 


For the longitudinal stratified conditional search probability in Figure 3.5, the selected 
search boxes also start with highest pi and slowly work to the left as pi lowers. This follows 
the slightly unexpected results from the discrete search of Chapter 2 possibly from the factor 
of 1 - qi being 93.7% in the western side and 37% in the east. This coupled with the pi 
being a factor of 10 higher in the east than the west leads it to be a higher probability of 
detection even though the sensor performs worse there. The probability of the target not 
being in the search area pQ = .1 ^ 0.1458. 
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Figure 3.5. Stratified Location and Longitudinal Stratified Search Results 
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Remaining probability of detection and boxes search based on stratified pi and 
linear stratified qi 


For the corner stratified conditional search probability in Figure 3.6, the selected boxes 
start in the southeast due to the juxtaposition of high pi and low qi. The path moves 
mostly along discretized radial arcs centered on the southeast corner. Which box is chosen 
via the Bayesian search algorithm is dependent on when 1 - qi will push a probability of 
detection higher for a lower p, . The probability of the target not being in the search area 
po = -1 ^ 0.1710. 


Figure 3.6. Stratified Location and Corner Stratified Search Density Results 



Remaining probability of detection and boxes search based on stratified pi and 
corner stratified q^ 


For the longitudinal stratified p, location probability density, the following results were 
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found based on running the seareh for up to 500 glimpses or a total probability of deteetion 
at 95%. A data set using a step size of 5% is plotted in Figure 3.7 along with associated fit 
lines as calculated from MATLAB’s fit algorithm based on 

f{x) = a 


Each of the fit lines is color coded to the data set with black for uniform, red for longitudinal 
stratification, and green for corner stratification. As the number of glimpses increases, the 
probability approaches Pa = -9, with the uniform case increasing faster than the other two. 
The associated parameters for the best fit lines are presented in Table 3.1. 

Figure 3.7. Stratified Location Cumulative Probability and Fit Lines 


F*d vs Number of Olimpses mid Exponential Fitted Curves 



Cumulative probabilities of detection based on the number of glimpses. Included 
are the fit best fit lines 


The results are very similar to results from Chapter 2 as the same decision algorithm is 
used; however, the exponential growth approaches Pa = -9 the initial probability that the 
target is in the search area A. 
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Table 3.1. Parameter Variation for p, for Longitudinal Stratification 


Parameter 

a 

b 


SSE 

RMSE 

P stratiQuniform 

P strut tQlong 

P strut t^corner 

0.8938 

0.90437 

0.89101 

-0.015155 

-0.0097729 

-0.0131 

0.99013 

0.99496 

0.97996 

0.0086678 

0.0053481 

0.017876 

0.0093465 

0.017237 

0.031514 


Best fit line parameters for a longitudinally stratified location probability density 
of the form f{x) = a - 


3.3.2 Gaussian Corner Location 

Figures 3.8, 3.9, and 3.10 show the remaining loeation probability densities for the three 
eonditional seareh probability densities after 50 glimpses. Also shown are the boxes that 
were seleeted to seareh by the Bayesian seareh algorithm. 


For the uniform eonditional seareh probability in Figure 3.8, unsurprisingly, the seleeted 
boxes started with highest p, and slowly works outward as p, shifts higher probability to 
other boxes. A seeond look begins oeeurring in those higher probability boxes prior to 
searehing all boxes. 


Figure 3.8. Gaussian Corner Location and Uniform Search Density Results 



Remaining probability of detection and boxes search based on Gaussian p,- 
and uniform qi 


For the longitudinal stratified eonditional seareh probability in Figure 3.9, the seleeted boxes 
also started with highest p,- and slowly work outward as p,- lowers. This follows the slightly 
unexpeeted result from Chapter 2 sinee there are better probabilities of deteetion in some 
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areas of lower target loeation probability. As the loeation probability shifts to other boxes, 
the higher pi inereases the likelihood of the target being located to the west of the initial 
location locus. From the Remainder Probability graph in Figure 3.9, the higher location 
probability remains dispersed in the northwest corner albeit non-uniformly. 


Figure 3.9. Gaussian Corner and Longitudinal Stratified Search Density Re¬ 
sults 



Remaining probability of detection and boxes search based on stratified pi and 
linear stratified qt 


For the corner stratified conditional search probability in Figure 3.10, the selected boxes 
start in the northeast due to the high pi overriding the high qi. The Remainder Probability 
graph in Figure3.10 shows some areas 2-3 times greater than other areas, especially around 
the original locus. Some boxes have been searched twice or three times before every box 
has been searched once. Since the location probability distribution pi remains identical to 
that from Chapter 2 (other than the being 90% of the original), the overall patterns of the 
search are identical to the results from Chapter 2. 
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Figure 3.10. Gaussian Corner and Corner Stratified Search Density Results 



Remaining probability of detection and boxes search based on stratified pi and 
corner stratified qi 


For the Gaussian corner pi location probability density, the following results were found 
based on running the search for up to 500 glimpses or a total probability of detection at 
95%. A data set using a step size of 5% is plotted in Figure 3.11 along with associated fit 
lines as calculated from MATLAB’s fit algorithm based on 

fix) = a - 


Each of the fit lines are color coded to the data set with black for uniform, red for longitudinal 
stratification, and green for corner stratification. As the number of glimpses increase, the 
probability approaches Pa = -9, with the uniform case increasing faster than the other two. 
The associated parameters for the best fit lines are presented in Table 3.2. 
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Figure 3.11. Gaussian Corner Cumulative Probability and Fit Lines 


F*d vs Number of Olimpses and Exponential Fitted Cnrves 



Cumulative probabilities of detection based on the number of glimpses. Included 


are the fit best fit lines 


Table 3.2. Parameter Variation for pi for Corner Gaussian Distribution 


Parameter 

a 

b 


SSE 

RMSE 

P GC iQ^uni f orm 
PGCiQlong 
PGC->Qcorner 

0.89231 

0.88774 

0.88725 

-0.03706 

-0.022588 

-0.01732 

0.98424 

0.97697 

0.9747 

0.0068884 

0.014581 

0.018924 

0.019563 

0.028462 

0.032424 


Best fit line parameters for a corner Gaussian location probability density of the 
form f{x) = a - 


The results are very similar to results from Chapter 2 as the same deeision algorithm is 
used. 


3.3.3 Cumulative Detection Probability 

The total probability of deteetion produeed by the Bayesian seareh algorithm after 50 
glimpses is shown in Table 3.3. The Gaussian p,’s have higher probabilities, whieh is not 
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unexpected since they have a much higher concentration on a few points. This allows for 
multiple searches over a smaller number of points instead of 50 different boxes. The other 
result seems to be that higher probability results when pi matches up to 1 - qi as it does 
when pstratified has a higher probability of detection when paired with qcomerjrat instead 
of qiongjrat- This follows the same pattern as the Discrete Search Method but lower overall 
since there is some probability that the target is outside of A. The increase to po after 50 
glimpses is shown in Table 3.4. 


Table 3.3. Pd after 50 glimpses 


Pi 

Uniform Long. Stratified 

Corner Stratified 

Uniform 

0.315 

0.36495 

0.40415 

Long. Stratified 

0.45818 

0.31435 

0.41507 

Corner Stratified 

0.51441 

0.39562 

0.36028 

Gaussian 

0.50339 

0.47669 

0.53271 

Gaussian Side 

0.64693 

0.50164 

0.58658 

Gaussian Corner 

0.71978 

0.57289 

0.49499 


Total probability of detection after 50 glimpses based on various combinations 


Table 3.4. po after 50 glimpses 


Pi \ 

Uniform Long. Stratified 

Corner Stratified 

Uniform 

0.1460 

0.1575 

0.1678 

Long. Stratified 

0.1846 

0.1458 

0.1710 

Corner Stratified 

0.2059 

0.1655 

0.1563 

Gaussian 

0.2014 

0.1911 

0.2140 

Gaussian Side 

0.2832 

0.2007 

0.2419 

Gaussian Corner 

0.3569 

0.2341 

0.1980 


Probability of target not in A after 50 glimpses based on various combinations 
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CHAPTER 4: 
Continuous Search 


4.1 Theory 

The following section outlines the basic concepts and theory behind a continuous exhaustive 
search. 

Terms 

A - Total search area (no searches conducted outside of A) 

Vs - Searcher’s velocity 
vt - Target’s velocity 

R - Range of the sensor. Cookie cutter sensor that will detect the target if the range 
is less than R. R = R(vs) 
w - Sweep width, w = 2R 

The sensor is a "cookie cutter" sensor that has 100% chance of detecting the target if the 
distance from the searcher to the target is less than R. In [6] both R and both vs were 
independent of each other and varied accordingly. Here, R is a function of vs which more 
closely resembles real-world sensor performance, especially in SONAR sensors. 

4.1.1 General Concepts 

The basic assumptions for this continuous search model are [6] [9] [13]: 

• Single target undergoing Brownian diffusion 

• Two searchers with "cookie-cutter" sensors 

• Exhaustive search using one of three search patterns 

• Sensors radius is velocity dependent R(vs) 

• Each searcher results independent of the other 

• No false positives 

• Non-adaptive search. 
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As shown in Figure 4.1 the sweep width for a moving seareher {w = 2R) will eover an 
area searehed A{t) = wvt. The exhaustive seareh eovers all of A without overhanging 
any previously searehed area or any area outside of A. The total time to coverage then is 



Figure 4.1. Search Area and Sweep Width. Source: [9]. 



Search Area (A) as it relates to Sweep Width {w) 


4.1.2 Model Overview 

Based on the work done by H. Wang and H. Zhou [9] and Z. Lukens’s thesis [6], the two 
searchers will move along one of three set paths: the inline ladder, the multipath ladder and 
the inline spiral search. For each path setup, 100,000 independent Monte Carlo simulations 
are conducted with the target uniformly distributed over A. The independent targets move 
randomly with a Brownian motion with a given diffusion rate. 

The searchers will move along the given path independently of each other without coordi¬ 
nation. An input parameter of the initial separation of the two searcher is given so that the 
searchers don’t overlap throughout the search. The separation in the inline searches was 
estimated at 4Lx where is half the width of A [6]. 

4.2 Problem and Model 

The model is thus two searchers (aircraft, surface vessels, subs, etc.) of similar capability 
performing the inline spiral in Figure 4.2, the inline ladder search in Figure 4.3, and the 
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multipath ladder search in Figure 4.4. 


Figure 4.2. Spiral Search. Source: [9]. 


Figure 4.3. Inline Ladder Search. Source: [9]. 




Figure 4.4. Multipath Ladder Search Pattern. Source: [9]. 



4.2.1 Parameter 

As an extension of previous work done [6], this chapter will explore the search patterns 
with a sensor that obeys Equation 4.1. 

RiVs) = (4.1) 

Rq is the performance of the sensor at = 0 and V is the exponential decay factor based 
on the sensor performance at speed. Each sensor, in reality, will behave differently based 
on its characteristics. Eor instance, a towed SONAR array requires a minimum speed in 
order to be drooped correctly and at the proper depth, whereas a hull-mounted SONAR 
does not have this minimum effect and is only limited at higher speed by self-noise. Eurther 
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compounding the real life sensors is that, at higher speeds, the self noise has non-linear 
effects based on cavitation above a threshold. To simplify the issue, the sensor is assumed 
to be operating in the range where Equation 4.1 holds. Each search velocity will be varied 
as Vs 6 (8,9,10,11,12) and the velocity decay factor from V e (5,6,7,8,9,10) giving 150 
different combinations for each search pattern. The initial search radius, Rq = 2.5 nm, will 
be set based on the largest search radius in [6]. 


4.3 Results 

In Tables 4.1, 4.2, and 4.3, the time in hours to achieve a probability of detection at 90% 
(Pd = 90%) are shown for V = 9 the inline ladder search, the multipath ladder, and the 
spiral search, respectively. 

Table 4.1. Time to Pd = 90% {V = 9) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

6.8662 

6.9782 

7.1117 

7.5250 

7.7083 

9 

7.5069 

7.8241 

7.7043 

7.8639 

7.9255 

Ship 2 10 

9.7362 

9.0589 

8.7202 

8.4741 

8.3054 

11 

10.1011 

10.0333 

9.2198 

8.9486 

8.6854 

12 

10.0425 

10.1511 

10.0968 

9.2826 

8.9569 


Inline ladder search 
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Table 4.2. Time to Pd = 9Q%{V = 9) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

8.1866 

8.2470 

8.2123 

8.1351 

8.0883 

9 

9.3043 

8.9871 

8.7202 

8.3385 

8.1969 

Ship 2 10 

10.3288 

9.5669 

9.2282 

8.7452 

8.4140 

11 

10.7790 

10.0333 

9.4231 

9.0164 

8.6854 

12 

10.4768 

10.5311 

9.7168 

9.1197 

8.7940 


Multipath ladder search 


Table 4.3. Time to Pd = 90%(y = 9) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

9.6390 

9.3043 

9.6515 

10.2367 

10.3139 

9 

9.8330 

9.4100 

9.3129 

9.7621 

9.8254 

Ship 2 10 

10.4135 

9.4822 

9.3129 

9.5587 

9.6625 

11 

10.5078 

10.0333 

9.5587 

9.6943 

9.7168 

12 

10.3682 

11.6710 

11.5625 

10.6397 

10.5854 


Spiral search 


In a marked difference from the results found in [6], faster search speed does not result in a 
quicker time to reach Pd = .9 in all cases. For instance, the lowest time to detection for the 
inline ladder search is with ai = 8 and V 2 = 8, for multipath ladder vi = 12 and V 2 = 8, and 
for spiral search vi = 9 and V 2 = S. In fact, keeping V 2 = S and varying vi over the range 
causes an increase and then decrease. Similar results happen on some other variations but 
as can be seen it is not consistent from one speed to another. Another noteworthy result is 
that varying the speed of ship 1 has less of an effect normally than varying ship 2’s speed. 
An example of this can be seen in the multipath ladder. Raising v\ from 8 to 12 knots keeps 
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the time to detection relatively constant around 8 hours, while raising vi from 8 to 12 knots 
increase the time to detect to approximately 8.5 hours. 

Similar to the results found in [6], the inline ladder search has the lowest time to detect, 
followed by the multipath ladder, and finally the spiral search for a given set of parameters. 
As expected, increasing V only lowers the time to detection. This follows [6] as increasing 
V raises the search radii without changing the speed. Direct comparison to other results 
from [6] is problematic. This work varies three of the four parameters from Equation 4.1, 
directly varying vi and V 2 and indirectly and In [6], and are varied directly. If 
Equation 4.1 is constructed with a V and Rq so that and 7?2 equal values from [6], then 
results will be the same as the search algorithm work in the same manner. 

A similar presentation for the various searches is outlined below in Eigures 4.5-4.10. 

Spiral Search 

Figure 4.5. Pd vs Time Figure 4.6. Pd vs Time 
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Multipath Ladder 

Figure 4.7. Pd vs Time Figure 4.8. Pd vs Time 



0 5 10 15 20 25 0 5 10 15 20 25 

Time (hrs) Time (hrs) 


Inline Ladder Search 



Figure 4.9. Pd vs Time 

Figure 4.10. Pd vs Time 
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As can be seen from the figures, there is an approach to approximately Pd = ^ - ■ An 

increase in V, as stated before increases P(vs) for a given Vs. This obviously hastens Pd to 
higher level more quickly. 
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CHAPTER 5: 

Conclusions and Future Work 


Based on the results from Chapters 2-4 the following sections summarize the conclusions drawn. 

5.1 Results and Conclusions 

5.1.1 Discrete and Bayesian Searches 

Both the Discrete and Bayesian searches follow the greedy algorithm and give similar 
results. In fact, if pq = 0, then the results will be identical to within rounding error. As 
expected, the more concentrated pi density, the higher the probability to detect Pj) and the 
lower number of glimpse needed to reach it. This even has the capacity to override higher 
conditional miss probability in that same concentrated area. Both search systems provide 
predictable, repeatable results, but the Bayesian search provides a more useful output. By 
modifying the pi at each step, not only can po > 0, but faulty conjectures that created the 
Pi be identified by scrutinizing the changing pi and increasing po. This provides a useful 
product by the algorithm at any point, the pi have been calculated if the number of glimpses 
have been exceeded and can be used at a later date. For a discrete search, the output from 
the algorithm is just the failure to detect on previous searches. 

5.1.2 Continuous Searches 

As Vs increases, Po has a non-intuitive result. It varies in a non-linear, way some times 
increasing or decreasing as Vs changes. Another non-intuitive result is that vi has a different 
change in Po than V 2 . The last result, an expected one, is that increasing V increases the Po 
by increasing R. 


5.2 Future Work 

5.2.1 Discrete and Bayesian Searches 

Improvements to the Discrete and Bayesian model in future work could include: 
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• Multiple Pi tested simultaneously updating all based on multiple conjectures 

• Multiple qi tested simultaneously updating all based on multiple exclusive sensors 

• Time dependent qi such as qi = as greedy algorithm not the optimal solution 

• Time dependent pi for a moving target. 


5.2.2 Continuous Searches 

Improvements to the continuous search in future work could include: 

• Different search patterns 

• Non-cookie cutter sensor with some non-uniform probability of detection 

• Game theory approach with a target counter-detecting and evading. 
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APPENDIX A: 
Figures 


A.l Discrete 


Uniform Location and Uniform Search Density Results 



Figure A.l. Remaining probability of detection and boxes search based on 
uniform pi and uniform qi 


Uniform Location and Longitudinal Stratified Search Density Results 



Figure A.2. Remaining probability of detection and boxes search based on 
uniform pi and linear stratified qi 
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Uniform Location and Corner Stratified Search Density Results 



Figure A.3. Remaining probability of detection and boxes search based on 
uniform pi and corner stratified qi 


Uniform Location Cumulative Probability and Fit Lines 


VS Number of Olimpses mid Exponential Fitted Curves 



Figure A.4. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Corner Stratified Location and Uniform Search Density Results 



Figure A.5. Remaining probability of detection and boxes search based on 
corner stratified pi and uniform qi 


Corner Stratified Location and Longitudinal Stratified Search Density Results 



Figure A.6. Remaining probability of detection and boxes search based on 
corner stratified pi and linear stratified qi 
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Corner Stratified Location and Corner Stratified Search Density Results 



Figure A.7. Remaining probability of detection and boxes search based on 
corner stratified pi and corner stratified qi 


Corner Stratified Location Cumulative Probability and Fit Lines 


F*d vs Number of Olimpses mid Exponential Fitted Curves 



Figure A.8. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Gaussian Location and Uniform Search Density Results 



Figure A.9. Remaining probability of detection and boxes search based on 
Gaussian pi and uniform qi 


Gaussian Location and Longitudinal Stratified Search Density Results 



Figure A.10. Remaining probability of detection and boxes search based on 
Gaussian pi and linear stratified qi 
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Gaussian Location and Corner Stratified Search Density Results 



Figure A.11. Remaining probability of detection and boxes search based on 
Gaussian pi and corner stratified qi 


Gaussian Location Cumulative Probability and Fit Lines 



Figure A.12. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Side Gaussian Location and Uniform Search Density Results 



Figure A.13. Remaining probability of detection and boxes search based on 
side Gaussian pi and uniform qi 


Side Gaussian Location and Longitudinal Stratified Search Density Results 



Long. 



Long. 


Figure A.14. Remaining probability of detection and boxes search based on 
side Gaussian pi and linear stratified qi 
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Side Gaussian Location and Corner Stratified Search Density Results 



Figure A.15. Remaining probability of detection and boxes search based on 
side Gaussian pi and corner stratified qi 


Side Gaussian Location Cumulative Probability and Fit Lines 



Figure A.16. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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A.2 Bayesian 

Uniform Location and Uniform Search Density Results 



Figure A.17. Remaining probability of detection and boxes search based on 
uniform pi and uniform qi 


Uniform Location and Longitudinal Stratified Search Density Results 



Figure A.18. Remaining probability of detection and boxes search based on 
uniform pi and linear stratified qi 
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Uniform Location and Corner Stratified Search Density Results 



Figure A.19. Remaining probability of detection and boxes search based on 
uniform pi and corner stratified qi 


Uniform Location Cumulative Probability and Fit Lines 


VS Number of Olimpses mid Exponential Fitted Curves 



Figure A.20. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Corner Stratified Location and Uniform Search Density Results 



Figure A.21. Remaining probability of detection and boxes search based on 
corner stratified pi and uniform qi 


Corner Stratified Location and Longitudinal Stratified Search Density Results 



Figure A.22. Remaining probability of detection and boxes search based on 
corner stratified pi and linear stratified qi 
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Corner Stratified Location and Corner Stratified Search Density Results 



Figure A.23. Remaining probability of detection and boxes search based on 
corner stratified pi and corner stratified qi 


Corner Stratified Location Cumulative Probability and Fit Lines 


F*d vs Number of Olimpses mid Exponential Fitted Curves 



Figure A.24. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Gaussian Location and Uniform Search Density Results 



Figure A.25. Remaining probability of detection and boxes search based on 
Gaussian pi and uniform qi 


Gaussian Location and Longitudinal Stratified Search Density Results 



Figure A.26. Remaining probability of detection and boxes search based on 
Gaussian pi and linear stratified qi 


53 






































































































































































































Gaussian Location and Corner Stratified Search Density Results 



Figure A.27. Remaining probability of detection and boxes search based on 
Gaussian pi and corner stratified qi 


Gaussian Location Cumulative Probability and Fit Lines 



Figure A.28. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Side Gaussian Location and Uniform Search Density Results 



Figure A.29. Remaining probability of detection and boxes search based on 
side Gaussian pi and uniform qi 


Side Gaussian Location and Longitudinal Stratified Search Density Results 



Figure A.30. Remaining probability of detection and boxes search based on 
side Gaussian pi and linear stratified qi 
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Side Gaussian Location and Corner Stratified Search Density Results 



Long. 



Long. 
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Figure A.31. Remaining probability of detection and boxes search based on 
side Gaussian pi and corner stratified qi 


Side Gaussian Location Cumulative Probability and Fit Lines 


vs Number of Olimpses mid Exponential Fitted Curves 



Figure A.32. Cumulative probabilities of detection based on the number of 
glimpses. Included are the fit best fit lines 
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Figure A.97. Pu vs Time Figure A.98. Pd vs Time 
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Figure A.99. Pd vs Time Figure A.100. Pd vs Time 
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Figure A. 124. Pd vs Time 
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Figure A.185. Pd vs Time Figure A.186. Pd vs Time 
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Figure A.187. Pd vs Time Figure A.188. Pd vs Time 
























































































































































































































































































Probability of Detection Probability of Detection 


Spiral 


A.3.3 Spiral Search 




Figure A.333. Pj) vs Time Figure A.334. Pd vs Time 
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Figure A.335. Pd vs Time 


Figure A.336. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.337. Pd vs Time Figure A.338. Pd vs Time 




Figure A.339. Pd vs Time 


Figure A.340. Pd vs Time 
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Probability of Detection Probability of Detection 
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Figure A.343. Pd vs Time 
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Figure A.342. Pd vs Time 
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Figure A.344. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.345. Pd vs Time 



Figure A.347. Pd vs Time 
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Figure A.348. Pd vs Time 


135 






Probability of Detection Probability of Detection 
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Figure A.349. Pd vs Time 


Figure A.350. Pd vs Time 



Figure A.351. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.353. Pd vs Time 



Figure A.355. Pd vs Time 



Figure A.354. Pd vs Time 
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Figure A.356. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.357. Pd vs Time Figure A.358. Pd vs Time 




Figure A.359. Pd vs Time 


Figure A.360. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.361. Pd vs Time Figure A.362. Pd vs Time 



Figure A.363. Pd vs Time 



Figure A.364. Pd vs Time 


139 
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Figure A.367. Po vs Time 
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Figure A.368. Pd vs Time 
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Probability of Detection Probability of Detection 
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Figure A.369. Pd vs Time 



Figure A.371. Pd vs Time 
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Figure A.372. Pd vs Time 
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Probability of Detection Probability of Detection 
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Figure A.373. Pd vs Time 


Figure A.374. Pd vs Time 



Figure A.375. Pd vs Time 
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Figure A.376. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.377. Pb vs Time Figure A.378. Pd vs Time 
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Figure A.379. Pd vs Time 


Figure A.380. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.381. Pd vs Time 



Figure A.383. Pd vs Time 


R,=0.920, R2=0.558, V^=8, V2=12, V=8 

1-^^-T-- =— 

0.9 - 
0.8 - 
I 0.7 - 

lo.6- 

□ 

0 0.5 - ■ 

lo.4- 

£l 

0 

[L 0.3 - /' 

0.2 - 
0.1 -/ 

ot- ^^^^- 

0 5 10 15 20 25 

Time (hrs) 

Figure A.382. Pd vs Time 



Figure A.384. Pd vs Time 
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Probability of Detection Probability of Detection 
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Figure A.385. Pd vs Time 
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Figure A.387. Pd vs Time 
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Figure A.386. Pd vs Time 



Time (hrs) 


Figure A.388. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.389. Pd vs Time 


Figure A.390. Pd vs Time 



Figure A.391. Pd vs Time 



Figure A.392. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.393. Pd vs Time 



Figure A.395. Pd vs Time 
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Figure A.396. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.397. Pd vs Time 



Figure A.399. Pd vs Time 



Figure A.398. Pd vs Time 
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Figure A.400. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.401. Pd vs Time Figure A.402. Pd vs Time 
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Figure A.404. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.405. Pd vs Time 



Figure A.407. Pd vs Time 
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Figure A.408. Pd vs Time 
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Figure A.409. Pd vs Time 
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Probability of Detection Probability of Detection 



Figure A.413. Pd vs Time 
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Figure A.416. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.417. Pd vs Time Figure A.418. Pd vs Time 




Figure A.419. Pd vs Time 


Figure A.420. Pd vs Time 
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Probability of Detection Probability of Detection 
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Figure A.424. Pd vs Time 
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Figure A.428. Pd vs Time 
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Figure A.431. Pd vs Time 
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Figure A.432. Pd vs Time 
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Figure A.434. Pd vs Time 



Figure A.435. Pd vs Time 


R,= 0 . 797 , R 2 = 0 . 519 , V^= 8 , V 2=1 1 , V =7 


0.9 


0.8 


I 0.7 


0 0.5 


lO'l 

£l 

0 

[L 0.3 


0.2 


0 '-^^^^-' 

0 5 10 15 20 25 

Time (hrs) 

Figure A.436. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.437. Pb vs Time Figure A.438. Pd vs Time 
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Figure A.440. Pd vs Time 
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Probability of Detection Probability of Detection 
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Figure A.444. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.445. Pd vs Time Figure A.446. Pd vs Time 




Figure A.447. Pd vs Time Figure A.448. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.449. Pd vs Time 


Figure A.450. Pd vs Time 



Figure A.451. Pd vs Time 



Figure A.452. Pd vs Time 


161 






Probability of Detection Probability of Detection 



Figure A.453. Pd vs Time 
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Figure A.456. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.457. Pd vs Time Figure A.458. Pd vs Time 




Figure A.459. Pd vs Time 


Figure A.460. Pd vs Time 


163 






Probability of Detection Probability of Detection 
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Probability of Detection Probability of Detection 



Figure A.465. Pd vs Time 
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Figure A.468. Pd vs Time 
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Figure A.470. Pd vs Time 


1 — 
0.9 - 
0.8 - 
0.7 - 
0.6 - 
0.5 - 
0.4 - 
0.3 - 
0.2 - 
0.1 


R^=0.920, R2=0.632, V^=8, V^sll, V=8 


0^ 

0 


10 15 20 25 


Time (hrs) 

Figure A.471. Pd vs Time 


R,=0.920, R2=0.558, V^=8, V2=12, V=8 


0.9 - 


1 0 . 6 - 
□ 

0 0.5 - 

l0-4- 

£l 

0 

[L 0.3 - 
0.2 - , 
0.1 7 


0^ 

0 


Time (hrs) 

Figure A.472. Pd vs Time 


166 






Probability of Detection Probability of Detection 



Figure A.473. Pd vs Time 
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Figure A.474. Pd vs Time 
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Probability of Detection Probability of Detection 




Figure A.477. Pd vs Time Figure A.478. Pd vs Time 
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Figure A.480. Pd vs Time 
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APPENDIX B: 
Tables 


B.l Discrete 

Table B.l. Parameter Variation for pi for Uniform Distribution 


Parameter 

b 

r2 

SSE 

RMSE 

P uni f ormiQuni f orm 

0.01061 

0.99242 

0.006608 

0.027096 

PuniformiQl ong 

0.009988 

0.99935 

0.00049844 

0.007442 

P uni f ornii^corner 

0.01198 

0.99795 

0.0015533 

0.013137 


Best fit line parameters for a uniformly distributed location probability density of 
the form f{x) - \ 


Table B.2. Parameter Variation for pi for Corner Stratification 


Parameter 

b 

r2 

SSE 

RMSE 

P corner str at iQuni f orm 

0.01321 

0.99908 

0.00065853 

0.008554 

P corner strat ^Qlong 

0.00896 

0.9992 

0.00066179 

0.0085751 

P corner strut Incomer 

0.01119 

0.99808 

0.0013515 

0.012254 


Best fit line parameters for a corner stratified distribution location probability 
density of the form f{x) = 1 - 
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Table B.3. Parameter Variation for pi for Centered Gaussian 


Parameter 

b 

r2 

SSE 

RMSE 

PGauss^^uniform 

0.0147 

0.99781 

0.0014554 

0.012717 

PGaussiQlong 

0.01315 

0.99609 

0.0026167 

0.017051 

PGauss^^corner 

0.01589 

0.9963 

0.00234 

0.016124 


Best fit line parameters for a centered Gaussian distribution location probability 
density of the form f{x) = 1 - 


Table B.4. Parameter Variation for pi for Side Gaussian 


Parameter 

b 

r2 

SSE 

RMSE 

PGauss^^uniform 

0.02405 

0.9944 

0.0027756 

0.017561 

PGausSfQlong 

0.01428 

0.99225 

0.0047039 

0.0228621 

PGauss jQcorner 

0.01895 

0.98981 

0.0054218 

0.024544 


Best fit line parameters for a offset Gaussian distribution location probability den¬ 
sity of the form f{x) = 1 - 


B.2 Bayesian 

Table B.5. Parameter Variation for pi for Uniform Distribution 


Parameter 

a 

b 

r2 

SSE 

RMSE 

P St rat uni form 

0.90643 

-0.01139 

0.99645 

0.004106 

0.015103 

Pstrat'>Qlong 

0.89218 

-0.011534 

0.98964 

0.0104 

0.024037 

P stratt^corner 

0.8876 

-0.014151 

0.99384 

0.0058881 

0.018086 


Best fit line parameters for a uniformly distributed location probability density of 
the form f{x) - a - 
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Table B.6. Parameter Variation for pi for Corner Stratification 


Parameter 

a 

b 

r2 

SSE 

RMSE 

P strut iQuni form 

0.88918 

-0.018268 

0.976 

0.0175 

0.031181 

Pstrut ■>^long 

0.89054 

-0.012323 

0.98204 

0.016665 

0.030427 

P strut corner 

0.89926 

-0.011051 

0.99147 

0.008804 

0.022116 


Best fit line parameters for a corner stratified distribution location probability 
density of the form f{x) = a - 


Table B.7. Parameter Variation for pi for Centered Gaussian 


Parameter 

a 

b 

r2 

SSE 

RMSE 

P strut iQunif orm 

0.88971 

-0.017384 

0.9847 

0.012037 

0.025859 

P strut t^long 

0.88274 

-0.016095 

0.98036 

0.01587 

0.029693 

P strut tQcorner 

0.88519 

-0.019317 

0.98268 

0.01265 

0.02651 


Best fit line parameters for a centered Gaussian distribution location probability 
density of the form f{x) = a - 


Table B.8. Parameter Variation for pi for Side Gaussian 


Parameter 

a 

b 

r2 

SSE 

RMSE 

P St r at uni form 

0.88957 

-0.029116 

0.97623 

0.012426 

0.026275 

P strut tQlong 

0.88243 

-0.01768 

0.96822 

0.023033 

0.035772 

P strut ■>Qcorner 

0.88301 

-0.023588 

0.96697 

0.019693 

0.0330771 


Best fit line parameters for a offset Gaussian distribution location probability den¬ 
sity of the form f{x) - a - 
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B.3 Continuous 

B.3.1 Inline Ladder Search 

Table B.9. Inline Ladder Search: Time to Pd = 90%(y = 5) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

14.8081 

14.5798 

15.6685 

16.1428 

16.6395 

9 

15.7111 

15.6258 

16.6558 

16.8813 

17.6232 

Ship 2 10 

18.5017 

18.0867 

17.8148 

18.2146 

18.7484 

11 

18.2434 

19.5095 

19.5958 

18.8285 

20.0086 

12 

18.1569 

19.1277 

20.5486 

21.9053 

21.3588 


Inline Ladder Search 


Table B.IO. Ini 

ine Ladder Search: 

Time to Pd 

= 90%(y 

= 6) 




Ship 1 



Velocities 

8 

9 

10 

11 

12 

8 

12.7567 

12.2523 

12.1515 

12.4612 

13.1644 

9 

13.3414 

12.6024 

12.3187 

12.5611 

13.1930 

Ship 2 10 

14.0466 

13.0154 

12.6810 

13.2800 

13.5936 

11 

13.7793 

15.2770 

14.1587 

14.0189 

14.6239 

12 

13.9657 

15.0961 

16.4268 

15.6398 

15.5540 


Inline Ladder Search 
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Table B.ll. Inline Ladder Search: Time to Pd = 9Q%{V = 7) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

10.9647 

10.6280 

10.4097 

10.2503 

10.2367 

9 

11.7028 

11.2848 

10.7237 

10.5538 

10.5661 

Ship 2 10 

12.7877 

11.5763 

11.0378 

10.8235 

10.7182 

11 

12.2397 

12.4420 

11.3968 

11.0259 

10.9716 

12 

11.8331 

12.0358 

12.8213 

11.4023 

11.4783 


Inline Ladder Search 


Table B.12. Inline Ladder Search: Time to Pd = 9Q%{V = 8) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

8.7756 

8.6460 

8.8498 

8.7901 

8.8294 

9 

10.1282 

9.7988 

9.2987 

9.2396 

9.0239 

Ship 2 10 

11.4149 

10.2606 

9.8759 

9.5892 

9.3740 

11 

11.1874 

11.3872 

10.2884 

9.8888 

9.6073 

12 

10.9687 

11.0076 

11.0854 

10.1519 

9.7240 


Inline Ladder Search 


175 





Table B.13. Inline Ladder Search: Time to Pd = 9Q%{V = 9) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

6.8662 

6.9782 

7.1117 

7.5250 

7.7083 

9 

7.5069 

7.8241 

7.7043 

7.8639 

7.9255 

Ship 2 10 

9.7362 

9.0589 

8.7202 

8.4741 

8.3054 

11 

10.1011 

10.0333 

9.2198 

8.9486 

8.6854 

12 

10.0425 

10.1511 

10.0968 

9.2826 

8.9569 


Inline Ladder Search 


Table B.14. Inline Ladder Search: Time to Pd = 9Q%{V = 10) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

6.6247 

6.4570 

6.3438 

6.4058 

6.5204 

9 

6.5861 

6.3278 

6.3438 

6.4924 

6.8038 

Ship 2 10 

7.2954 

6.9782 

7.0840 

7.0118 

7.1582 

11 

9.1759 

8.6565 

8.2237 

7.8774 

7.6543 

12 

9.2135 

9.2844 

9.0009 

8.3631 

8.1504 


Inline Ladder Search 
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B.3.2 Multipath Ladder 


B.3.3 Multipath Ladder Search 

Table B.15. Multipath Ladder Search: Time to Pd = 90%(y = 5) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

14.6171 

14.9854 

15.6398 

16.1907 

16.5752 

9 

14.9427 

15.8392 

16.3410 

16.9964 

17.6940 

Ship 2 10 

18.2155 

17.4714 

17.5716 

18.2338 

18.7098 

11 

18.2913 

19.0874 

18.5983 

18.9148 

19.9443 

12 

18.3240 

19.1792 

20.6129 

20.6580 

21.1080 


Multipath Ladder Search 


Table B.16. Multipath Ladder Search: Time to Pd = 9Q%{V = 6) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

12.7024 

12.2912 

12.2072 

13.0803 

13.1644 

9 

12.9135 

12.4857 

12.4580 

12.9206 

13.2073 

Ship 2 10 

13.7122 

12.7089 

12.7089 

13.2600 

13.7510 

11 

14.0389 

14.0988 

13.6395 

14.2585 

14.5953 

12 

13.9371 

15.6255 

15.8545 

15.2106 

15.5254 


Multipath Ladder Search 
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Table B.17. Multipath Ladder Search: Time to Pd = 9Q%{V = 7) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

11.4414 

10.9266 

10.4994 

10.2503 

10.3888 

9 

11.8222 

11.1654 

10.7686 

10.6212 

10.7942 

Ship 2 10 

12.3390 

11.3070 

10.9481 

10.8235 

11.0729 

11 

12.4083 

11.7002 

11.1270 

10.8910 

11.3263 

12 

12.0611 

12.3145 

11.5037 

11.0729 

11.630 


Multipath Ladder Search 


Table B.18. Multipath Ladder Search: Time to Pd = 9Q%{V = 8) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

10.1501 

9.7165 

9.2987 

9.0398 

8.9850 

9 

10.7870 

10.2105 

9.7476 

9.2895 

9.1406 

Ship 2 10 

11.2226 

10.5813 

9.9400 

9.5392 

9.45180 

11 

11.5370 

10.7379 

10.0886 

9.7890 

9.6073 

12 

11.0854 

11.2021 

10.3464 

9.8407 

9.6462 


Multipath Ladder Search 
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Table B.19. Multipath Ladder Search: Time to Pd = 9Q%{V = 9) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

8.1866 

8.2470 

8.2123 

8.1351 

8.0883 

9 

9.3043 

8.9871 

8.7202 

8.3385 

8.1969 

Ship 2 10 

10.3288 

9.5669 

9.2282 

8.7452 

8.4140 

11 

10.7790 

10.0333 

9.4231 

9.0164 

8.6854 

12 

10.4768 

10.5311 

9.7168 

9.1197 

8.7940 


Multipath Ladder Search 


Table B.20. Multipath Ladder Search: Time to Pd = 9Q%{V = 10) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

6.7825 

6.8444 

6.7668 

7.0983 

7.0873 

9 

6.9735 

7.4901 

7.4011 

7.5312 

7.3708 

Ship 2 10 

8.7756 

8.4585 

8.2470 

7.8774 

7.6543 

11 

9.9550 

9.1759 

8.6565 

8.3102 

8.0087 

12 

9.8514 

9.8514 

9.0009 

8.5757 

8.221 


Multipath Ladder Search 
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B.3.4 


Spiral 


B.3.5 Spiral Search 

Table B.21. Spiral Search: Time to Pd = 90%(y = 5) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

19.0117 

19.4895 

19.5176 

19.7109 

20.2015 

9 

21.6669 

20.4074 

20.8340 

20.8907 

20.0729 

Ship 2 10 

22.2220 

23.1092 

21.5638 

21.5909 

21.67388 

11 

21.8115 

22.8090 

25.0000 

22.5980 

22.65114 

12 

21.7960 

22.4003 

23.4676 

25.0000 

25.0000 


Spiral Search 


Table B.22. Spiral Search: Time to Pd = 90%{V = 6) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

13.8424 

14.3527 

14.3253 

14.5182 

16.0405 

9 

16.0252 

15.7918 

16.2205 

16.5751 

15.0818 

Ship 2 8 

18.3108 

17.6140 

17.1124 

17.3739 

17.5859 

11 

18.2326 

18.7518 

18.6519 

18.2126 

18.1725 

12 

18.3013 

18.9309 

19.5462 

19.8467 

18.9309 


Spiral Search 
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Table B.23. Spiral Search: Time to Pd = 90%{V = 7) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

11.5209 

11.8819 

12.1147 

12.1723 

12.5932 

9 

11.7625 

12.0013 

12.3839 

12.1723 

11.4023 

Ship 2 10 

15.1658 

13.7749 

13.5505 

13.7233 

14.0122 

11 

15.0720 

15.8476 

15.0383 

14.8697 

15.1017 

12 

14.8737 

16.1153 

16.5460 

15.9886 

15.7099 


Spiral Search 


Table B.24. Spiral Search: Time to Pd = 90%{V = 8) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

10.3616 

10.2929 

10.9019 

11.0875 

11.2021 

9 

10.4576 

10.2929 

10.6454 

10.6879 

10.5409 

Ship 2 10 

11.8638 

10.5813 

10.6454 

10.7878 

10.6575 

11 

12.2362 

13.2850 

11.9365 

11.8866 

11.7855 

12 

12.0189 

13.6137 

13.9637 

13.1469 

12.9524 


Spiral Search 
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Table B.25. Spiral Search: Time to Pd = 90%{V = 9) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

9.6390 

9.3043 

9.6515 

10.2367 

10.3139 

9 

9.8330 

9.4100 

9.3129 

9.7621 

9.8254 

Ship 2 10 

10.4135 

9.4822 

9.3129 

9.5587 

9.6625 

11 

10.5078 

10.0333 

9.5587 

9.6943 

9.71684 

12 

10.3682 

11.6710 

11.5625 

10.6397 

10.5854 


Spiral Search 


Table B.26. Spiral Search: Time to Pd = 90%{V = 10) 


Velocities 

8 

9 

Ship 1 

10 

11 

12 

8 

9.3062 

8.6524 

8.6699 

9.5221 

9.4262 

9 

9.2981 

8.7815 

8.5642 

8.8296 

9.2135 

Ship 2 10 

9.6215 

8.8814 

8.5642 

8.5699 

8.8592 

11 

9.6087 

9.1759 

8.5699 

8.6565 

8.7883 

12 

9.3553 

9.7805 

9.0009 

8.7174 

8.788 


Spiral Search 
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APPENDIX C: 
Codes 


C.l 

Discrete 

C.1.1 

Main Code 

1 ele 


2 elear 



3 discrete_search_prob 

4 discrete_loc_probs2 

5 close all 

6 %p-uniform 

7 dis_pUni_qUni 

8 dis_pUni_qlongstrat 

9 dis_pUni_qbothstrat 

10 %p-long stratified 

11 dis_plongstrat_qUni 

12 dis_plongstrat_qlongstrat 

13 dis_plongstrat_qbothstrat 

14 %p-both stratified 

15 dis_pbothstrat_qUni 

16 dis_pbothstrat_qlongstrat 

17 dis_pbothstrat_qbothstrat 

18 %p-normal eenter 

19 dis_pnormwide_qUni 

20 dis_normwide_qlongstrat 

21 dis_normwide_qbothstrat 

22 %p-normal side 

23 dis_pnormside_qUni 

24 dis_normside_qlongstrat 

25 dis_normside_qbothstrat 
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26 %p-normal corner 

27 dis_pnormc_qUni 

28 dis_normc_qlongstrat 

29 dis_normc_qbothstrat 

30 m 

31 

32 pd = [ Pd_uni_uni ( end ) , Pd_uni_ls ( end ) , Pd_uni_bs ( end ) ,Pd_ls_uni( 

end ) ,.. . 

33 Pd_ls_ls ( end ) , Pd_ls_bs ( end ) , Pd_bs_uni ( end ) , Pd_bs_ls ( end ) 

,Pd_bs_bs( end )... 

34 , Pd_nw_uni ( end ) , Pd_nw_ls ( end ) , Pd_nw_bs ( end ) , Pd_ns_uni ( 

end ) ,. . . 

35 Pd_ns_ls ( end ) , Pd_ns_bs ( end ) , Pd_ne_uni ( end ) , Pd_ne_ls ( end ) 

? • • • 

36 Pd_ne_bs ( end ) ] 

37 

38 m 

39 % elose all 

40 set(0,’defaulttextinterpreter’ , ’latex’) 

41 clear g 

42 [ ~ , Pd_uni_uni , k] = diserete_seareh2 ( P_uni , q_uni ,500 , 1) 

43 [ ~ , Pd_uni_ls ,k] = diserete_seareh2 ( P_uni , q_long_strat 

,500, 1 ); 

44 [ ~ , Pd_uni_bs ,k] = discrete_search2 ( P_uni , q_both_strat 

,500, 1 ); 

45 figure 

46 

47 stepl =10:round(.05*length(Pd_uni_uni)) : length(Pd_uni_uni) ; 

48 step2 = 10:round(.05*length(Pd_uni_ls)) : length(Pd_uni_ls) ; 

49 step3 = 10:round(.05*length(Pd_uni_bs)) : length(Pd_uni_bs) ; 

50 plot(stepl ,Pd_uni_uni(stepl ) , ’xk’ ,step2 ,Pd_uni_ls(step2) , ’xr 

’ , step3 , Pd_uni_bs (step3) , ’xg’) 
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51 g = fittype (’ 1-exp (b*x)dependentProbindependent’x 

52 

53 [ fit_bb 1 , gof_bb l]=fit(stepl ’ ,Pd_uni_uni(stepl) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

54 [ fit_bb2 , gof_bb2 ]=fit(step2 ’ ,Pd_uni_ls(step2) ’ ,g , ’StartPoint 

’ ,[0] , ’Upper’ ,[0]) ; 

55 [ fit_bb3 , gof_bb3 ] = fit (step3 ’ , Pd_uni_bs (step3) ’,g, ’StartPoint 

’ ,[0] , ’Upper’ ,[0]) ; 

56 hold on 

57 plot ( fit_bb 1 , ’k ’ ) 

58 plot ( fit_bb2 , ’ r ’ ) 

59 plot ( fit_bb3 , ’g ’ ) 

60 xlabel(’n’) 

61 ylabel ( ’$P_{D}$ ’ ) 

62 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

63 legend (’ (p_{ uniform } ,q_{ uniform })’,’(p_{ uniform } ,q_{ s tr at a 

}) ’ ,’ (p_{ uniform } ,q_{ strat b }) ’ ... 

64 ,sprintf(’l-e^{%.5f*x}’ ,fit_bbl .b) ,sprintf(’l-e^{%.5f*x} 

’ ,fit_bb2 .b) ,... 

65 sprintf(’l-e^{%.5f*x}’,fit_bb3.b),’looation’,’best’) 

66 print ( ’ Bayes_Uni_Fit ’ , ’-dpng ’ ) 

67 

68 [ ~ , Pd_ls_uni , k] = discrete_search2 ( P_long_strat , q_uni 

,500, 1 ); 

69 [ ~ , Pd_ls_ls ,k] = disorete_searoh2 ( P_long_strat , 

q_long_strat ,500 , 1 ); 

70 [ ~ , Pd_ls_bs ,k] = disorete_searoh2 ( P_long_strat , 

q_both_strat ,500, 1 ); 

71 figure 

72 stepl =10:round(.05*length(Pd_ls_uni)) : length(Pd_ls_uni) ; 

73 step2 = 10:round(.05*length(Pd_ls_ls)) : length(Pd_ls_ls) ; 
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74 step3 = 10:round(.05*length(Pd_ls_bs)) : length(Pd_ls_bs) ; 

75 plot(stepl ,Pd_ls_uni(stepl) ,’xk’ ,step2 ,Pd_ls_ls(step2) ,’xr’ , 

step3 ,Pd_ls_bs(step3) ,’xg’) 

76 [fit_bb4, gof_bb4 ] = fit (stepl ’ ,Pd_ls_uni(stepl ) ’ ,g, ’ StartPoint 

’ ,[0] , ’Upper’ ,[0]) ; 

77 [ fit_bb5 , gof_bb5 ]=fit(step2 ’ ,Pd_ls_ls(step2) ’,g, ’StartPoint’ 

,[0] , ’Upper’ ,[0]) ; 

78 [ fit_bb6 , gof_bb6 ]=fit(step3 ’ ,Pd_ls_bs(step3) ’,g, ’StartPoint ’ 

,[0] , ’Upper’ ,[0]) ; 

79 hold on 

80 plot ( fit_bb4 , ’k ’ ) 

81 plot ( fit_bb5 , ’ r ’ ) 

82 plot ( fit_bb6 , ’g ’ ) 

83 xlabel(’n’) 

84 ylabel ( ’$P_{D}$ ’ ) 

85 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

86 legend (’ (p_{ s tr at a } ,q_{ uniform })’,’ (p_{ s tr at a},q_{strat a 

})’ , ’(p_{strat a} ,q_{strat b })’ . .. 

87 ,sprintf(’l-e^{%.5f*x}’ ,fit_bb4.b) ,sprintf(’l-e^{%.5f*x} 

’ ,fit_bb5 .b) 

88 sprintf(’l-e^{%.5f*x}’ ,fit_bb6.b) ,’looation’ ,’best’) 

89 print ( ’ Bayes_Long_Fit ’ , ’-dpng ’ ) 

90 

91 [ ~ , Pd_bs_uni , k] = diserete_seareh2 ( P_both_strat , q_uni 

,500, 1 ); 

92 [ ~ , Pd_bs_ls ,k] = disorete_searoh2 ( P_both_strat , 

q_long_strat ,500 , 1 ); 

93 [ ~ , Pd_bs_bs ,k] = disorete_searoh2 ( P_both_strat , 

q_both_strat ,500, 1 ); 

94 figure 

95 stepl =10:round(.05*length(Pd_bs_uni)) : length(Pd_bs_uni) ; 

96 step2 = 10:round(.05*length(Pd_bs_ls)) : length(Pd_bs_ls) ; 
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97 step3 = 10:round(.05*length( Pd_bs_bs )): length ( Pd_bs_bs ) ; 

98 plot(stepl, Pd_bs_uni (stepl) ,’xk’ ,step 2 , Pd_bs_ls (step2) ,’xr’ , 

step 3 , Pd_bs_bs (step3),’xg’) 

99 [ fit_bb7 , gof_bb7 ]=fit(stepl Pd_bs_uni (stepl ) ’ ,g, ’ StartPoint 

’ ,[0] , ’Upper’ ,[0]) ; 

100 [ fit_bb8 , gof_bb8 ] = fit (step2 ’ ,Pd_bs_ls(step2) ’ ,g, ’ StartPoint ’ 

,[0] , ’Upper’ ,[0]) ; 

101 [ fit_bb9 , gof_bb9 ] = fit (step3 ’ , Pd_bs_bs (step3) ’ ,g, ’ StartPoint ’ 

,[0] , ’Upper’ ,[0]) ; 

102 hold on 

103 plot(fit_bb7 , ’k’) 

104 plot ( fit_bb8 , ’ r ’ ) 

105 plot ( fit_bb9 , ’g ’ ) 

106 Xlabel (’ n ’ ) 

107 ylabel ( ’$P_{D}$ ’ ) 

108 1111 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

109 legend (’ (p_{ s tr at b } ,q_{ uniform })’,’ (p_{ s tr at b},q_{strat a 

})’,’(P_{strat b},q_{strat b})’... 

110 ,sprintf(’l-e^{%.5f*x}’ ,fit_bb7 .b) ,sprintf(’l-e^{%.5f*x} 

’ ,fit_bb8 .b) 

111 sprintf(’l-e^{%.5f*x}’ ,fit_bb9.b) ,’looation’ ,’best’) 

112 print ( ’ Bayes_Both_Fit ’ , ’-dpng ’ ) 

113 

114 [ ~ ,Pd_nw_uni , k] = discrete_search2 (P_norm_wide , q_uni 

,500, 1 ); 

115 [ ~ , Pd_nw_ls ,k] = disorete_searoh2 (P_norm_wide , 

q_long_strat ,500 , 1 ); 

116 [ ~ ,Pd_nw_bs ,k] = disorete_searoh2 (P_norm_wide , 

q_both_strat ,500, 1 ); 

117 figure 

118 stepl =10: round(.05*length( Pd_nw_uni )):length( Pd_nw_uni) ; 

119 step2 = 10:round(.05*length( Pd_nw_ls )):length( Pd_nw_ls ) ; 
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120 step3 = 10: round(. 05* length (Pd_nw_bs)) : length (Pd_nw_bs) ; 

121 plot ( step 1 , Pd_nw_uni (stepl) ,’xk’ ,step 2 , Pd_nw_ls (step2) ,’xr’ , 

step 3 , Pd_nw_bs (step3),’xg’) 

122 [ fit_bb 1 0 , gof_bb 10 ] = f i t ( step 1 ’ , Pd_nw_uni (stepl) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

123 [ fit_bb 1 1 , gof_bb 11 ] = f i t ( step 2 ’ , Pd_nw_ls (step2) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

124 [ fit_bb 1 2 , gof_bb 12]=fit(step3 ’ , Pd_nw_bs (step3) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

125 hold on 

126 plot ( fit_bb 1 0 , ’k’) 

127 plot ( fit_bb 11 ,’r’) 

128 plot ( fit_bb 1 2 , ’g’) 

129 Xlabel (’ n ’ ) 

130 ylabel ( ’$P_{D}$ ’ ) 

131 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

132 legend (’ (p_{norm } ,q_{ uniform (p_{norm } ,q_{ strat a})’,’( 

p_ { norm } , q_ { s t r a t b }) ’ . . . 

133 ,sprintf(’l-e^{%.5f*x}’ ,fit_bbl0.b) ,sprintf(’l-e^{%.5f*x 

} ’ ,fit_bbll .b) 

134 sprintf(’l-e^{%.5f*x}’ ,fit_bbl2.b) ,’looation’ ,’best’) 

135 print ( ’ Bayes_Norm_Fit ’ , ’-dpng ’ ) 

136 

137 [ ~ , Pd_ns_uni , k] = discrete_search2 ( P_norm_side , q_uni 

,500, 1 ); 

138 [ ~ , Pd_ns_ls ,k] = disorete_searoh2 ( P_norm_side , 

q_long_strat ,500 , 1 ); 

139 [ ~ , Pd_ns_bs ,k] = disorete_searoh2 ( P_norm_side , 

q_both_strat ,500, 1 ); 

140 figure 

141 stepl =10:round(.05*length(Pd_ns_uni)) : length(Pd_ns_uni) ; 

142 step2 = 10:round(.05*length(Pd_ns_ls)) : length(Pd_ns_ls) ; 


188 


143 step3 = 10:round(.05*length( Pd_ns_bs )):length( Pd_ns_bs ) ; 

144 plot ( step 1 , Pd_ns_uni (stepl) ,’xk’ ,step 2 ,Pd_ns_ls(step2) ,’xr’ , 

step 3 , Pd_ns_bs (step3),’xg’) 

145 [ fit_bb 1 3 , gof_bbl3 ] = f it ( step 1 ’ , Pd_ns_uni (stepl) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

146 [ fit_bb 1 4 , gof_bb 14 ] = f i t ( step 2 ’ ,Pd_ns_ls(step2) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

147 [ fit_bb 1 5 , gof_bb 15]=fit(step3 ’ , Pd_ns_bs (step3) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

148 hold on 

149 plot ( fit_bb 1 3 , ’k’) 

150 plot(fit_bbl4,’r’) 

151 plot ( fit_bb 15 , ’g ’ ) 

152 Xlabel (’ n ’ ) 

153 ylabel ( ’$P_{D}$ ’ ) 

154 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

155 legend (’ (p_{ norm b }, q_{ uniform })’,’(p_{ norm b},q_{strat a})’ 

,’(p_{norm b},q_{strat b})’... 

156 ,sprintf(’l-e^{%.5f*x}’ ,fit_bbl3 .b) ,sprintf(’l-e^{%.5f*x 

} ’ ,fit_bbl4 .b) 

157 sprintf(’l-e^{%.5f*x}’ ,fit_bbl5 .b) ,’looation’ ,’best’) 

158 print ( ’ Bayes_NormSide_Fit ’ , ’-dpng ’ ) 

159 

160 [ ~ , Pd_no_uni , k] = diserete_seareh2 ( P_norm_oorner , q_uni 

,500, 1 ); 

161 [ ~ , Pd_no_ls ,k] = disorete_searoh2 ( P_norm_oorner , 

q_long_strat ,500 , 1 ); 

162 [ ~ , Pd_no_bs ,k] = disorete_searoh2 ( P_norm_oorner , 

q_both_strat ,500, 1 ); 

163 figure 

164 stepl =10:round(.05*length(Pd_ne_uni)) : length(Pd_ne_uni) ; 

165 step2 = 10:round(.05*length(Pd_ne_ls)) : length(Pd_ne_ls) ; 


189 


166 step3 = 10:round(.05*length( Pd_nc_bs )):length( Pd_nc_bs ) 

167 plot ( step 1 , Pd_nc_uni (stepl) ,’xk’ ,step 2 ,Pd_ne_ls(step2) ,’xr’ , 

step 3 , Pd_ne_bs (step3),’xg’) 

168 [ fit_bb 1 6 , gof_bbl6 ] = f it ( step 1 ’ , Pd_ne_uni (stepl) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

169 [ fit_bb 17 , gof_bb 17 ] = f i t ( step 2 ’ ,Pd_ne_ls(step2) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

170 [ fit_bb 1 8 , gof_bb 18]=fit(step3 ’ , Pd_ne_bs (step3) ’,g, ’ 

StartPoint ’ ,[0] , ’ Upper ’ , [0]) ; 

171 hold on 

172 plot ( fit_bb 1 6 , ’k’) 

173 plot(fit_bbl7,’r’) 

174 plot ( fit_bb 18 , ’g ’ ) 

175 Xlabel (’ n ’ ) 

176 ylabel ( ’$P_{D}$ ’ ) 

177 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

178 legend (’ (p_{ norm e }, q_{ uniform })’,’(p_{ norm o},q_{strat a})’ 

,’(p_{norm e},q_{strat b})’... 

179 ,sprintf(’l-e^{%.5f*x}’ ,fit_bbl6.b) ,sprintf(’l-e^{%.5f*x 

} ’ ,fit_bbl7 .b) 

180 sprintf(’l-e^{%.5f*x}’,fit_bbl8.b),’looation’,’best’) 

181 print ( ’ Bayes_NormCorner_Fit ’ , ’-dpng ’ ) 

182 format shortG 

183 b = [fit_bbl .b, fit_bb2 .b, fit_bb3 .b, fit_bb4 .b, fit_bb5 .b, fit_bb6 

. b , fit_bb7 . b , . . . 

184 fit_bb8 .b, fit_bb9 .b, fit_bbl0 .b, fit_bbl 1 .b, fit_bbl2 .b, 

fit_bbl3 .b, fit_bbl4 .b ,. . . 

185 fit_bbl5 .b,fit_bbl6 .b,fit_bbl7 .b,fit_bbl8 .b] 

186 r2 =[gof_bb 1 . rsquare , gof_bb2 . rsquare , gof_bb3 . rsquare , gof_bb4 . 

rsquare , . . . 

187 gof_bb5 . rsquare , gof_bb6 . rsquare , gof_bb7 . rsquare , gof_bb8 . 

rsquare , . . . 
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gof_bb9 . rsquare , gof_bbl0 . rsquare , gof_bbl 1 . rsquare , 
gof_bbl2 . rsquare , . . . 

gof_bbl3 . rsquare , gof_bbl4 . rsquare , gof_bbl5 . rsquare , 
gof_bbl6 . rsquare , . . . 
gof_bbl7 . rsquare , gof_bbl8 . rsquare ] 

C.1.2 Search Probabilities qi 

format short 
format eompaet 
Area = 1000; 
x=-5:5;y=x; 

[X,Y] = meshgrid(x,y) ; 

m 

q0 = .3; 

q_uni=q0*ones(10,10); 
figure 

subplot (1 ,3,1) 

imagese ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , q_uni , [0 ,1 ]) 

set(gea, ’YDir’, ’normal’) 

ob_q=colorbar ; 

title (eb_q , ’q_{i}’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e ( ’ Uniform ’ ) 
xlabel ( ’Long ’ ) 
ylabel ( ’ Lat ’ ) 

m 

% 

q0 = .07; 

qx = 1:10; 

qy = ones (10,1) ; 
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q_long_strat=qy*qx*q0*.9; 
subplot (1 ,3,2) 

image sc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] ,q_long_strat ,[0 ,max( 
q_long_strat (:))]) 
set(gca, ’YDir’, ’normal’) 
cb_q=colorbar ; 
title (cb_q , ’q_{i}’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title(’Strat. a’) 
xlabel ( ’Long ’ ) 
ylabel ( ’ Lat ’ ) 

m 

q0 = .01; 
qx = (l:10) ; 
qy=(l:10) ’; 

q_both_strat=qy*qx*q0*.9; 
subplot (1 ,3 ,3) 

imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] ,q_both_strat ,[0 ,max( 
q_both_strat (:))]) 
set(gca, ’YDir’, ’normal’) 
cb_q=colorbar ; 
title (cb_q , ’q_{i}’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title(’Strat. b’) 
xlabel ( ’Long ’ ) 
ylabel ( ’ Lat ’ ) 

m 

print ( ’ Discrete_Miss_Probablities ’ , ’-dpng ’ ) 
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C.1.3 Location Probabilities pi 

% clear 

close all 
clc 

format short 
format compact 
Area = 1000; 

%create the grid positions 

p0 = 0; 

x= -5:5 ;y=x; 

[X,Y] = meshgrid(x,y) ; 

m 

P_uni = ones(10,10) ; 

P_uni=(l-p0)*P_uni /sum( P_uni (:) ) ; 
figure 

subplot (2,3,1) 

imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , P_uni , [0 ,1 ]) 

set(gca, ’YDir’, ’normal’) 

cb_p=colorbar ; 

title (cb_p , ’p_{i}’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title ( ’ Uniform ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat ’ ) 

m 

Px = 1:10; 

Py=ones (10,1) ; 

P_long_strat=Py*Px; 

P_long_strat=(l-pO)*P_long_strat /sum( P_long_strat (:) ) ; 
subplot (2,3,2) 

imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] ,P_long_strat ,[0 ,max( 
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P_long_strat (:))]) 
set(gca, ’YDir’, ’normal’) 
cb_p=colorbar ; 
title (cb_p , ’p_{i } ’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title ( ’ Stratified a ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat ’ ) 

m 

Px = 1:10; 

Py = (l:10) ’; 

P_both_strat=Py*Px; 

P_both_strat=(l-pO)*P_both_strat /sum( P_both_strat (:) ) ; 
% figure 
subplot (2 ,3 ,3) 

imagese ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] ,P_both_strat ,[0 ,max( 
P_both_strat (:))]) 
set(gca, ’YDir’, ’normal’) 
ob_p=eolorbar ; 
title (cb_p , ’p_{i}’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11e ( ’ Stratif i e d b’) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat ’ ) 

m 

Px = normpdf ([-4.5 :1:4.5] ,0 ,2) ; 

Py = normpdf([-4.5:1:4.5],0,5) ’; 

P_norm_wide=Py*Px; 

P_norm_wide=(l -pO) *P_norm_wide / sum( P_norm_wide (:) ) ; 
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subplot (2,3,4) 

imagesc([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , P_norm_wide , [ 0 , max( P_norm_wide 

(:))]) 

set(gca, ’YDir’, ’normal’) 
cb_p=colorbar ; 
title (cb_p , ’p_{i } ’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e ( ’ Normal a ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat ’ ) 

m 

Px = normpdf ([-4.5 :1:4.5] ,3 ,2) ; 

Py = normpdf([-4.5:1:4.5],3,2) ’; 

P_norm_corner=PyH<Px; 

P_norm_corner=(l -pO) *P_norm_corner /sum( P_norm_corner (:) ) ; 
subplot (2,3,5) 

image sc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , P_norm_corner , [0 , max( 
P_norm_corner (:))]) 
set(gca, ’YDir’, ’normal’) 
cb_p=colorbar ; 
title (cb_p , ’p_{i}’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title ( ’ Normal b ’ ) 
xlabel ( ’Long. ’ ) 
ylabel( ’Lat’) 

m 

Px = normpdf ([-4.5 :1:4.5] ,3 ,2) ; 

Py = normpdf([-4.5:1:4.5],0,2) ’; 

P_norm_side=Py*Px; 
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P_norm_side=(l -pO) *P_norm_side / sum( P_norm_side (:) ) ; 
subplot (2,3,6) 

imagesc([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , P_norm_side , [ 0 , max( P_norm_side 

(:))]) 

set(gca, ’YDir’, ’normal’) 
cb_p=colorbar ; 
title (cb_p , ’p_{i}’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e ( ’ Normal c ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat ’ ) 

m 

print ( ’ Diserete_Loeation_Probablities ’ , ’-dpng ’ ) 

C.1.4 Sample code for the Uniform/Uniform case 

[ pe , n , imax , Pd_uni_uni ,k] = diserete_seareh2 ( P_uni , q_uni 
,50, .95 ); 
figure 

subplot (2,2,1) 

imagese ([ - 4.5,4.5],[-4.5,4.5],po,[0 ,max(po (:))]) 

set(goa, ’YDir’, ’normal’) 

ob_p=colorbar ; 

title (eb_p , ’p_{i}’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e ( ’ Remainder Prob . ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 

subplot (2,2,2) 
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imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] ,n,[0 , max(n (:))]) 

set(gca, ’YDir’, ’normal’) 

cb_p=colorbar ; 

title (cb_p , ’N’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e (’ Actual Search’) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 
subplot (2,2,3) 

imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , P_uni , [0 , max( P_uni (:) ) ]) 

set(gca, ’YDir’, ’normal’) 

cb_p=colorbar ; 

title (cb_p , ’p_{i}’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e ( ’ Uniform ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 


subplot (2,2,4) 

imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , q_uni , [0 , max( q_uni (:))]) 

set(gca, ’YDir’, ’normal’) 

cb_q=colorbar ; 

title (cb_q , ’q_{i}’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title ( ’ Uniform ’ ) 
xlabel ( ’Long. ’ ) 
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ylabel ( ’ Lat. ’ ) 

print ( ’Discrete_uni_uni ’ , ’-dpng’) 

C.1.5 Discrete Search Function 

function [ pd , n , imax , Pd_cum2 ,k] = diserete_search2 (p , q , nT , 
stop ) 

^Conducts a Discrete Search base on inputs . 

%p= Location probability density 
^cq=Conditional miss probabilit density 
%nT=max number of glimpses 

%stop = Probability of Detection to stop algorithm 
%Outputs 

%pd= updated location probability density 
%n=search matrix 
%imax = Last box searched 
%Pd_cum= total probabiliy of detection 
%k=number of glimpses 
if -exist ( ’ stop ’ ) 
stop=sum(p (:) ) ; 
elseif stop >sum(p (:) ) 
stop=sum(p (:) ) ; 

end 

p0=l-sum(p (:) ) ; 

n= zeros ( size (p)) ;%i n i tia 1 i ze look counter 
Pd_cum = 0; 

Pd_cum2=Pd_cum; 
pc=p; 
pc_0=p0; 
k=l; 

while k<=nT && Pd_cum2 ( end )< stop 
pd=pc .*(1 - q) .*q.^n; 

[pmax , imax (k) ] = max(pd (:) ) ; 
pmax; 
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[irow , icol] = ind2sub(size(pd) ,imax(k)) 

30 n(irow , icol )=n(irow , icol)+ 1; 

31 Pd_cum2(k+l)=Pd_cum2(k)+pmax; 

32 k=k+l; 

33 end 

34 pd=pc .*(1-q) .*q .^n ; 

35 end 


C.2 Bayesian 

C.2.1 Main Code 

1 clc 

2 clear all 

3 bayes_search_prob 

4 bayes_loc_probs2 

5 close all 

6 %p-uniform 

7 bayes_pUni_qUni 

8 p0_uni_uni = pc_0 

9 bayes_pUni_qlongstrat 

10 p0_uni_ls=pc_0 

11 bayes_pUni_qbothstrat 

12 p0_uni_bs = pc_0 

13 %p-long stratified 

14 bayes_plongstrat_qUni 

15 p0_ls_uni = pc_0 

16 bayes_plongstrat_qlongstrat 

17 p0_ls_ls =pc_0 

18 bayes_plongstrat_qbothstrat 

19 p0_ls_bs =pc_0 

20 %p-both stratified 

21 bayes_pbothstrat_qUni 
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pO_bs_uni=pc_0 

23 bayes_pbothstrat_qlongstrat 

24 pO_bs_ls =pc_0 

25 bayes_pbothstrat_qbothstrat 

26 p0_bs_bs=pc_0 

27 %p-normal center 

28 bayes_pnormwide_qUni 

29 p0_nw_uni = pc_0 

30 bayes_normwide_qlongstrat 

31 p0_nw_ls=pe_0 

32 bayes_normwide_qbothstrat 

33 p0_nw_bs=pc_0 

34 %p-normal side 

35 bayes_pnormside_qUni 

36 p0_ns_uni = pc_0 

37 bayes_normside_qlongstrat 

38 pO_ns_ls =pc_0 

39 bayes_normside_qbothstrat 

40 p0_ns_bs=pe_0 

41 %p-normal corner 

42 bayes_pnormc_qUni 

43 p0_nc_uni = pe_0 

44 bayes_norme_qlongstrat 

45 p0_nc_uni = pc_0 

46 bayes_normc_qbothstrat 

47 p0_nc_bs=pc_0 

48 

49 pd = [ Pd_uni_uni ( end ) , Pd_uni_ls ( end ) , Pd_uni_bs ( end ) ,Pd_ls_uni( 

end ) ,. . . 

50 Pd_ls_ls ( end ) , Pd_ls_bs ( end ) , Pd_bs_uni ( end ) , Pd_bs_ls ( end ) 

,Pd_bs_bs( end )... 

51 , Pd_nw_uni ( end ) , Pd_nw_ls ( end ) , Pd_nw_bs ( end ) , Pd_ns_uni ( 

end ) ,... 
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Pd_ns_ls( end ) ,Pd_ns_bs( end ) , Pd_nc_uni( end ) , Pd_nc_ls( end ) 

? • • • 

53 Pd_nc_bs ( end ) ] 

54 m 

55 

56 elose all 

57 set(0,’defaulttextinterpreter’ , ’latex’) 

58 [ ~ ,Pd_uni_uni ,k,~ ,~] = bayes_seareh ( P_uni , q_uni ,500 , 

.95 ); 

59 [ ~ , Pd_uni_ls ,k,~ ,~] = bayes_search ( P_uni , q_long_strat 

,500, .95 ); 

60 [ ~ ,Pd_uni_bs ,k,~ ,~] = bayes_search ( P_uni , q_both_strat 

,500, .95 ); 

61 figure 

62 stepl =10:round(.05*length(Pd_uni_uni)) : length(Pd_uni_uni) ; 

63 step2 = 10:round(.05*length(Pd_uni_ls)) : length(Pd_uni_ls) ; 

64 step3 = 10:round(.05*length(Pd_uni_bs)) : length(Pd_uni_bs) ; 

65 plot(stepl ,Pd_uni_uni(stepl ) , ’xk’ , step 2 ,Pd_uni_ls(step2) ,’xr 

’ , step3 , Pd_uni_bs (step3) , ’xg’) 

66 g = fittype (’ a-exp (b*x) ’,’ dependent ’,’ Prob ’,’ independent ’, ’x 

67 

68 [fit_bbl , gof_bb 1] = fit (stepl ’ ,Pd_uni_uni(stepl ) ’ ,g, ’ 

StartPoint ’ ,[0 ,0] , ’ Upper ’ , [ 1 ,0]) ; 

69 [ fit_bb2 , gof_bb2 ]=fit(step2 ’ ,Pd_uni_ls(step2) ’,g, ’StartPoint 

’ ,[0,0] , ’Upper’ ,[1 ,0]) ; 

70 [ fit_bb3 , gof_bb3 ]=fit(step3 ’ , Pd_uni_bs (step3) ’,g, ’StartPoint 

’ ,[0,0] , ’Upper’ ,[1 ,0]) ; 

71 hold on 

72 plot ( fit_bb 1 , ’k ’ ) 

73 plot ( fit_bb2 , ’ r ’ ) 

74 plot ( fit_bb3 , ’g ’ ) 

75 xlabel(’n’) 
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ylabel( ’$P_{D}$’) 
ylim([0,1]) 

t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 
Curves ’ ) 

legend( ’(p_{ uniform },q_ {uniform})’ ,’(p_ (uniform ),q_{strat a 
’ (P_{ uniform ),q_{strat b})’... 
,sprintf(’%.5f-e^{%.5f*x)’ ,fit_bbl .a,fit_bbl .b),sprintf( 
’%.5f-e^{%.5f*x) ’ ,fit_bb2 .a, fit_bb2 .b) 
sprintf ( ’%.5f-e^{%.5f*x) ’ ,fit_bb3 .a, fit_bb3 .b) , ’ loeation 
’ , ’best’ ) 

print ( ’ Bayes_Uni_Fit ’ , ’-dpng ’ ) 

m 

[ ~ , Pd_ls_uni ,k,~ ,~] = bayes_seareh ( P_long_strat , q_uni 

,500, .95 ); 

[ ~ , Pd_ls_ls ,k,~ ,~] = bayes_seareh ( P_long_strat , 

q_long_strat ,500, .95 ); 

[ ~ , Pd_ls_bs ,k,~ ,~] = bayes_seareh ( P_long_strat , 

q_both_strat ,500, .95 ); 
figure 

stepl =10:round(.05*length(Pd_ls_uni)) : length(Pd_ls_uni) ; 
step2 = 10:round(.05*length(Pd_ls_ls)) : length(Pd_ls_ls) ; 
step3 = 10:round(.05*length(Pd_ls_bs)) : length(Pd_ls_bs) ; 
plot(stepl ,Pd_ls_uni(stepl) ,’xk’ ,step2 ,Pd_ls_ls(step2) ,’xr’ , 
step3 ,Pd_ls_bs(step3) ,’xg’) 

[ fit_bb4 , gof_bb4 ] = fit (stepl ’ ,Pd_ls_uni(stepl ) ’ ,g, ’ StartPoint 
’ ,[1 ,0] , ’Upper’ ,[1 ,0]) ; 

[ fit_bb5 , gof_bb5 ]=fit(step2 ’ ,Pd_ls_ls(step2) ’,g, ’StartPoint’ 
,[1 ,0] , ’Upper’ ,[1 ,0]) ; 

[ fit_bb6 , gof_bb6 ]=fit(step3 ’ ,Pd_ls_bs(step3) ’,g, ’StartPoint ’ 
,[1 ,0] , ’Upper’ ,[1 ,0]) ; 
hold on 
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98 plot ( fit_bb4 , ’k ’ ) 

99 plot ( fit_bb5 , ’ r ’ ) 

100 plot ( fit_bb6 , ’g ’ ) 

101 Xlabel (’ n ’ ) 

102 ylabel ( ’$P_{D}$ ’ ) 

103 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

104 legend (’ (p_{ strat a } ,q_{ uniform (p_{ s tr at a},q_{strat a 

})’,’(P_{strat a},q_{strat b})’... 

105 ,sprintf(’%.5f-e^{%.5f*x}’ ,fit_bb4.a,fit_bb4.b),sprintf( 

’%.5f-e^{%.5f*x} ’ ,fit_bb5 .a, fit_bb5 .b) . 

106 sprintf(’%.5f-e^{%.5f*x}’,fit_bb6.a,fit_bb6.b),’looation 

’ , ’best’ ) 

107 print ( ’ Bayes_Long_Fit ’ , ’-dpng ’ ) 

108 

109 [ ~ , Pd_bs_uni , k ,~ ~] = bayes_seareh ( P_both_strat , q_uni 

,500, .95 ); 

110 [ ~ , Pd_bs_ls ,k,~ ,~] = bayes_seareh ( P_both_strat , 

q_long_strat ,500, .95 ); 

111 [ ~ ,Pd_bs_bs ,k,~ ,~] = bayes_seareh ( P_both_strat , 

q_both_strat ,500, .95 ); 

112 figure 

113 stepl =10:round(.05*length(Pd_bs_uni)) : length(Pd_bs_uni) ; 

114 step2 = 10:round(.05*length(Pd_bs_ls)) : length(Pd_bs_ls) ; 

115 step3 = 10:round(.05*length( Pd_bs_bs )):length( Pd_bs_bs ) ; 

116 plot ( step 1 , Pd_bs_uni (stepl ) , ’xk’ ,step 2 , Pd_bs_ls (step2) ,’xr’ , 

step 3 , Pd_bs_bs (step3),’xg’) 

117 [ fit_bb7 , gof_bb7 ]=fit(stepl ’, Pd_bs_uni (stepl ) ’ ,g, ’ StartPoint 

’ ,[1 ,0] , ’Upper’ ,[1 ,0]) ; 

118 [ fit_bb8 , gof_bb8 ] = fit (step2 ’ ,Pd_bs_ls(step2) ’ ,g, ’ StartPoint ’ 

,[1 ,0] , ’Upper’ ,[1 ,0]) ; 

119 [ fit_bb9 , gof_bb9 ]=fit(step3 ’ , Pd_bs_bs (step3) ’ ,g, ’ StartPoint ’ 

,[1 ,0] , ’Upper’ ,[1 ,0]) ; 
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120 


hold on 

121 plot ( fit_bb7 , ’k ’ ) 

122 plot ( fit_bb8 , ’ r ’ ) 

123 plot ( fit_bb9 , ’g ’ ) 

124 Xlabel (’ n ’ ) 

125 ylabel ( ’$P_{D}$ ’ ) 

126 1111 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

127 legend (’ (p_{ s tr at b } ,q_{ uniform (p_{ s tr at b},q_{strat a 

})’,’(P_{strat b},q_{strat b})’... 

128 ,sprintf(’%.5f-e^{%.5f*x}’,fit_bb7.a,fit_bb7.b),sprintf( 

’%.5f-e^{%.5f*x} ’ ,fit_bb8 .a,fit_bb8 .b) 

129 sprintf(’%.5f-e^{%.5f*x}’,fit_bb9.a,fit_bb9.b),’loeation 

’ , ’best’ ) 

130 print ( ’ Bayes_Both_Fit ’ , ’-dpng ’ ) 

131 

132 [ ~ , Pd_nw_uni , k ,~ ~] = bayes_seareh (P_norm_wide , q_uni 

,500, .95 ); 

133 [ ~ ,Pd_nw_ls ,k,~ ,~] = bayes_seareh (P_norm_wide , 

q_long_strat ,500, .95 ); 

134 [ ~ ,Pd_nw_bs ,k,~ ,~] = bayes_seareh (P_norm_wide , 

q_both_strat ,500, .95 ); 

135 figure 

136 stepl=10:round(.05*length( Pd_nw_uni )):length( Pd_nw_uni) ; 

137 step2 = 10:round(.05*length( Pd_nw_ls )):length( Pd_nw_ls ) ; 

138 step3 = 10:round(.05*length( Pd_nw_bs )):length( Pd_nw_bs) ; 

139 plot ( step 1 , Pd_nw_uni (stepl ) , ’xk’ ,step 2 , Pd_nw_ls (step2) ,’xr’ , 

step 3 , Pd_nw_bs (step3) ,’xg’) 

140 [ fit_bb 1 0 , gof_bb 10 ] = f i t ( step 1 ’ , Pd_nw_uni (stepl )’,g,’ 

StartPoint ’ ,[1 ,0] , ’ Upper ’,[1,0]); 

141 [ fit_bb 11 , gof_bb ll]=fit(step2 ’ , Pd_nw_ls (step2) ’,g, ’ 

StartPoint ’ ,[1 ,0] , ’ Upper ’ , [ 1 ,0]) ; 

142 [ fit_bb 1 2 , gof_bb 12]=fit(step3 ’ , Pd_nw_bs (step3) ’,g, ’ 
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StartPoint ’ ,[1 ,0] , ’ Upper ’ , [ 1 ,0]) ; 

143 hold on 

144 plot ( fit_bb 1 0 , ’k’) 

145 plot ( fit_bb 11 ,’r’) 

146 plot ( fit_bb 1 2 , ’g’) 

147 Xlabel (’ n ’ ) 

148 ylabel ( ’$P_{D}$ ’ ) 

149 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

150 legend (’ (p_{norm } ,q_{ uniform (p_{norm } ,q_{ strat a})’,’( 

p_ { norm } , q_ { s t r a t b }) ’ . . . 

151 ,sprintf(’%.5f-e^{%.5f*x}’ ,fit_bbl0.a,fit_bbl0.b) , 

sprintf(’%.5f-e^{%.5f*x}’ ,fit_bbll .a,fit_bbll .b) 

152 sprintf(’%.5f-e^{%.5f*x}’,fit_bbl2.a,fit_bbl2.b),’ 

loeation ’ , ’ best ’ ) 

153 print ( ’ Bayes_Norm_Fit ’ , ’-dpng ’ ) 

154 

155 [ ~ , Pd_ns_uni , k ,~ ~] = bayes_seareh ( P_norm_side , q_uni 

,500, .95 ); 

156 [ ~ , Pd_ns_ls ,k,~ ,~] = bayes_seareh ( P_norm_side , 

q_long_strat ,500, .95 ); 

157 [ ~ ,Pd_ns_bs ,k,~ ,~] = bayes_seareh ( P_norm_side , 

q_both_strat ,500, .95 ); 

158 figure 

159 stepl =10:round(.05*length(Pd_ns_uni)) : length(Pd_ns_uni) ; 

160 step2 = 10:round(.05*length(Pd_ns_ls)) : length(Pd_ns_ls) ; 

161 step3 = 10:round(.05*length( Pd_ns_bs )):length( Pd_ns_bs ) ; 

162 plot ( step 1 , Pd_ns_uni (stepl) ,’xk’ ,step 2 ,Pd_ns_ls(step2) ,’xr’ , 

step 3 , Pd_ns_bs (step3),’xg’) 

163 [ fit_bb 1 3 , gof_bb 13]=fit(stepl ’, Pd_ns_uni (stepl) ’,g, ’ 

StartPoint ’ ,[1 ,0] , ’ Upper ’,[1,0]); 

164 [ fit_bb 1 4 , gof_bb 14 ] = f i t ( step 2 ’ ,Pd_ns_ls(step2) ’,g, ’ 

StartPoint ’ ,[1 ,0] , ’ Upper ’,[1,0]); 
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165 [ fit_bb 1 5 , gof_bb 15] = fit ( step 3 ’ , Pd_ns_bs (step3) ’,g, ’ 

StartPoint ’ ,[1 ,0] , ’ Upper ’,[1,0]); 

166 hold on 

167 plot ( fit_bb 1 3 , ’k’) 

168 plot(fit_bbl4,’r’) 

169 plot ( fit_bb 1 5 , ’g’) 

170 Xlabel (’ n ’ ) 

171 ylabel ( ’$P_{D}$ ’ ) 

172 t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 

Curves ’ ) 

173 legend (’( p_{ norm b }, q_{ uniform })’,’(p_{ norm b},q_{strat a}) 

,’(p_{norm b},q_{strat b])’... 

174 ,sprintf(’%.5f-e^{%.5f*x}’ ,fit_bbl3 .a,fit_bbl3 .b) , 

sprintf(’%.5f-e^{%.5f*x}’ ,fit_bbl4.a,fit_bbl4.b) ,... 

175 sprintf(’%.5f-e^{%.5f*x}’,fit_bbl5.a,fit_bbl5.b),’ 

location ’ , ’ best ’ ) 

176 print ( ’ Bayes_NormSide_Fit ’ , ’-dpng ’ ) 

177 

178 [ ~ , Pd_nc_uni , k ,~ ,~] = bayes_search ( P_norm_corner , q_uni 

,500, .95 ); 

179 [ ~ , Pd_nc_ls , k ,~ ,~] = bayes_search ( P_norm_corner , 

q_long_strat ,500, .95 ); 

180 [ ~ ,Pd_nc_bs ,k,~ ,~] = bayes_search ( P_norm_corner , 

q_both_strat ,500, .95 ); 

181 figure 

182 stepl =10:round(.05*length(Pd_nc_uni)) : length(Pd_nc_uni) ; 

183 step2 = 10:round(.05*length(Pd_nc_ls)) : length(Pd_nc_ls) ; 

184 step3 = 10:round(.05*length( Pd_nc_bs )):length( Pd_nc_bs ) ; 

185 plot ( step 1 , Pd_nc_uni (stepl) ,’xk’ ,step 2 ,Pd_nc_ls(step2) ,’xr’ 

step 3 , Pd_nc_bs (step3) ,’xg’) 

186 [ fit_bb 1 6 , gof_bb 16]=fit(stepl ’, Pd_nc_uni (stepl) ’,g, ’ 

StartPoint ’ ,[1 ,0] , ’ Upper ’ , [ 1 ,0]) ; 

187 [ fit_bb 17 , gof_bbl7 ] = fit (step 2 ’ ,Pd_nc_ls(step2) ’,g, ’ 
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StartPoint ’ ,[1 ,0] , ’ Upper ’ , [ 1 ,0]) ; 

[ fit_bb 1 8 , gof_bb 18]=fit(step3 ’ , Pd_nc_bs (step3) ’,g, ’ 
StartPoint ’ ,[1 ,0] , ’ Upper ’ , [ 1 ,0]) ; 
hold on 

plot (fit_bbl6 , ’k’ ) 
plot ( fit_bb 17 , ’ r ’ ) 
plot (fit_bbl8 , ’g’) 
xlabel(’n’) 
ylabel( ’$P_{D}$’) 

t i 11 e ( ’ $P_{D} $ vs Number of Glimpses and Exponential Fitted 
Curves ’ ) 

legend (’ (p_{ norm e }, q_{ uniform })’,’(p_{ norm o},q_{strat a})’ 
,’(p_{norm e},q_{strat b})’... 

, sprintf ( ’%.5f-e''{%.5f*x} ’ ,fit_bbl6 .a, fit_bbl6 .b) , 
sprintf ( ’%.5f-e^{%.5f*x} ’ ,fit_bbl7 .a, fit_bbl7 .b) 
sprintf ( ’%.5f-e^{%.5f*x} ’ ,fit_bbl8 .a, fit_bbl8 .b) , ’ 
loeation ’ , ’ best ’ ) 

print ( ’ Bayes_NormCorner_Fit ’ , ’-dpng ’ ) 
format shortG 

b = [fit_bbl .b, fit_bb2 .b, fit_bb3 .b, fit_bb4 .b, fit_bb5 .b, fit_bb6 
. b , fit_bb7 . b , . . . 

fit_bb8 .b, fit_bb9 .b, fit_bbl0 .b, fit_bbl 1 .b, fit_bbl2 .b, 
fit_bbl3 .b, fit_bbl4 .b ,. . . 
fit_bbl5 .b, fit_bbl6 .b, fit_bbl7 .b, fit_bbl8 .b] 
r2 =[gof_bb 1 . rsquare , gof_bb2 . rsquare , gof_bb3 . rsquare , gof_bb4 . 
rsquare , . . . 

gof_bb5 . rsquare , gof_bb6 . rsquare , gof_bb7 . rsquare , gof_bb8 . 
rsquare , . . . 

gof_bb9 .rsquare , gof_bbl0 .rsquare , gof_bbl 1 . rsquare , 
gof_bbl2 . rsquare , . . . 

gof_bbl3 . rsquare , gof_bbl4 . rsquare , gof_bbl5 . rsquare , 
gof_bbl6 . rsquare , . . . 
gof_bbl7 . rsquare , gof_bbl8 . rsquare ] 
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[ ~ ,Pd_uni_bs ,k,~ ,~] = bayes_search ( P_uni , q_both_strat 
,500, 1); 

[ ~ , Pd_uni_bs_r ,k] = bayes_search_random ( P_uni , 

q_both_strat ,500, 1 ); 
figure 

plot (1: length(Pd_uni_bs) , Pd_uni_bs ,’xr’ ,l:length(Pd_uni_bs_r 
) , Pd_uni_bs_r, ’r ’ ) 

Wc 

for i = 1:9 

p0 = . 1 * i ; 

P_uni = ones (10,10) ; 

P_uni=(l - pO)*P_uni /sum( P_uni (:) ) ; 

[ ~ ,Pd_uni_bs_pO(i ,:) ,k,~ ,~] = bayes_search ( P_uni , 
q_both_strat ,500 , 1); 

end 

figure 
for i = 1:9 
hold on 

plot (1:25:500 , Pd_uni_bs_p0 (i ,1:25:500)) 

end 


C.2.2 Sample code for the Uniform/Uniform case 

[ pe , n , imax , Pd_uni_uni ,k , pe_0 , k_eheek ] = bayes_seareh ( P_uni 
, q_uni ,50 , .95 ) ; 
figure 

subplot (2,2,1) 

imagesc ([ - 4.5,4.5],[-4.5,4.5], pc,[0 , max (pc (:))]) 

set(gca, ’YDir’, ’normal’) 

cb_p=colorbar ; 

title (cb_p , ’p_{i}’ ) 

hold on 
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plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e (’ Remainder Prob . ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 


subplot (2,2,2) 

image sc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] ,n,[0 , max(n (:))]) 

set(gca, ’YDir’, ’normal’) 

cb_p=colorbar ; 

title (cb_p , ’N’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
t i 11 e (’ Actual Search’) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 
subplot (2,2,3) 

image sc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , P_uni , [0 , max( P_uni (:))]) 

set(gca, ’YDir’, ’normal’) 

cb_p=colorbar ; 

title (cb_p , ’p_{i}’ ) 

hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title ( ’ Uniform ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 


subplot (2,2 ,4) 

imagesc ([ - 4.5 ,4.5] ,[ - 4.5 ,4.5] , q_uni , [0 , max( q_uni (:))]) 
set(gca, ’YDir’, ’normal’) 
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cb_q=colorbar ; 
title (cb_q , ’q_{i}’ ) 
hold on 

plot(X,Y, ’k’) ; 
plot(Y,X, ’k’) ; 
title ( ’ Uniform ’ ) 
xlabel ( ’Long. ’ ) 
ylabel ( ’ Lat . ’ ) 

print( ’Bayes_uni_uni ’ , ’-dpng ’ ) 

C.2.3 Bayes Search Function 

funetion [ pd , n , imax , Pd_eum ,k , pe_0 , k_eheck ] = bayes_seareh ( 
p,q,nT, stop ) 

%Conduots a Bayesian Seareh base on inputs . 

%p=Loeation probability density 
%q=Conditional miss probabilit density 
%nT=max number of glimpses 

%stop = Probability of Deteetion to stop algorithm 
%Outputs 

%pd= up dated loeation probability density 

%n=seareh matrix 

%imax = Last box searehed 

%k=number of glimpses 

%pe_0 =final probability tgt is not in seareh area 
%k_eheek=game theory deteetion glimpse 
if -exist ( ’ stop ’ ) 
stop=sum(p (:) ) ; 
elseif stop >sum(p (:) ) 
stop=sum(p (:) ) ; 

end 

eheek =0; 

n= zeros ( size (p)) ;%i n i tia 1 i ze look eounter 
Pd_oum = 0; 
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pc_0=l-sum(p (:) ) ; 

pc=p; 

k=l; 

tgt_pos = r and sample (1:101 ,1 ,true ,[p(:) ; pc_0 ]) ; 

[ tgt_row , tgt_col]=ind2sub(size(p) ,tgt_pos) ; 
while k<=nT && Pd_cum( end )< stop 
pd=pe .*(1 -q) ; 

[pmax , imax (k) ] = max(pd (:) ) ; 

[irow , ieol] = ind2sub(size(pd) ,imax(k)) ; 
if eheek==0 

if irow==tgt_row && ieo 1 == tgt_eo 1 

eheek = binornd(l,l -q(irow , ieol )) ; 
if check==l 

k_eheck=k; 

end 

end 

end 

pe(irow , ieol )=pe(irow , ieol ) .*q(irow , ieol) ; 

po=po/(l -pmax) ; 

po_0=pc_0/(l -pmax) ; 

n(irow , ieol )=n(irow , ieol) + 1; 

Pd_oum(k+l)=Pd_oum(k) +(1 -Pd_cum(k) )*pmax ; 

k=k+ 1; 

end 

if cheek==0 

k_eheck=NaN ; 

end 

pd=pc; 
end 


C.3 Continuous 

Code is as in [6] with the exeeption that Ri = and R 2 = Roe~'^^^^ 
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